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FOREWORD 

This  report  to  the  Office  of  Naval  Research  describes  both  theoreti¬ 
cal  and  experimental  efforts  carried  out  under  contract  N00014-80-C-0498  dur¬ 
ing  the  period  June  6,  1980  through  May  31,  1981.  The  program  was  carried  out 
at  the  Rockwell  International  Science  Center,  and  was  managed  by  Dr.  Paul  R. 
Newman.  The  principal  investigators  were  Mr.  Mark  D.  Ewbank  and  Dr.  Newman. 
Prof.  Walter  A.  Harrison  of  Stanford  University  was  a  consultant  on  the 
theoretical  aspects  of  tne  program.  Valuable  contributions  were  also  made  by 
Dr.  Pochi  Veil,  Mr.  Randolph  L.  Hall  and  Dr.  M.  Khoshnevisan  of  the  Science 
Center.  The  contract  monitor  for  the  Office  of  Naval  Research  was 
Dr.  George  Weight. 
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1.0  INTRODUCTION 


1.1  Program  Objectives 

This  program  has  several  long  range  objectives.  The  first  and  perhaps 
foremost  is  to  achieve  a  fuller  understanding  at  the  microscopic  level,  of  the 
physics  leading  to  the  electro-optic  (E-0)  effect  in  solids.  A  further  goal  is 
to  generate  a  completely  generalized  predictive  methodology  which  takes  as  its 
input  atomic  structure  and  elemental  composition  and  successfully  predicts  the 
electro-optic  response  for  the  material. 

Additionally,  efforts  in  electro-optic  device  development  programs  a^e 
currently  reaching  performance  levels  limited  not  by  device  design,  but  rather 
by  materials  characteri sties.  It  is  hoped  that  some  of  the  future  experimental 
efforts  of  this  program  will  lead  to  the  successful  identification,  growth,  and 
characterization  of  new  high  performance  electro-optic  materials,  so  badly 
needed  by  this  emerging  technology. 

1 . 2  Previous  Work 

Work  that  is  being  carried  out  under  this  program  has  evolved  from  and 
is  a  natural  consequence  of  a  Rockwell  International  IR&D  program  (Project  864 
"Optical  Materials")  at  the  Science  Center.  During  FY  1979  and  1980,  work  was 
begun  in  conjunction  with  Prof.  W.A.  Harrison  of  Stanford  University,  on  the 
application  of  his  "Bond-Orbital  Model or  chemical -bond  approach  to  the 
prediction  of  the  optical  dielectric  susceptibility  of  solids;  specifically  the 
indices  of  refraction  of  solids.  After  several  initial  successes  in  predicting 
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the  first-order  susceptibility,  efforts  evolved  towards  second-order  "perturba¬ 
tion-response"  functions^  such  as  the  strain-optic  effect  and  the  dc  electro¬ 
optic  (Pockels)  effect.  That  portion  of  the  effort  which  is  specifically 
oriented  toward  the  electro-optic  effect  is  the  basis  for  this  program  which  was 
officially  funded  on  June  1,  1980.  The  remainder  of  the  work  on  the  strain- 
optic  effect,  the  anomalous  dispersion  in  the  birefringence,  and  other  basic 
optical  properties  of  solids  will  continue  to  be  carried  out  at  the  Science 
Center  under  IRSO  funding.  Clearly,  there  will  be  many  synergistic  advantages 
both  theoretically  and  experimentally  in  having  such  a  comprehensive  approach  to 
the  generic  topic  of  the  interaction  between  light  and  solids. 


1 . 3  Accompl i shments 

There  have  been  several  significant  achievements  during  the  first  phase 
of  this  project.  The  theoretical  work  has  progressed  to  first  separating  the 
physical  sources  of  the  electro-optic  effect  into  those  which  result  primarily 
from  electronic  effects  and  those  associated  with  lattice  dynamics  or  "ionic- 
displacive  effects."  The  "Bond-Orbital  Model"  was  then  modified  so  that  the 
electronic  contribution  to  the  electro-optic  effect  could  be  calculated  simply 
and  straightforwardly  using  "universal  atomic"  parameters  and  sums  over  chemical 
bonds  in  the  crystallographic  unit  cell.  This  theory  was  tested  against  Te02 
and  predicted  an  identically  zero  electronic  contribution  to  the  electro-optic 
tensor  in  agreement  with  Kleinman's^  symmetry. 

The  ionic  displacive  part  was  next  to  be  investigated.  Here,  a  further 
distinction  was  made  between  relative  motions  between  ions  which  do  not  result 
in  a  change  in  the  unit  cell  dimensions,  and  those  which  do.  The  approach. 
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although  not  completely  evolved,  has  been  to  seek  a  mechanism  for  calculating 
the  microscopic  atomic  " spring  constants"  and  then  relating  these  to  some  more 
easily  measured  bulk  properties  such  as  the  Reststrahl  frequency  or  bulk  elastic 
constants.  It  is  hoped  that  this  will  enable  predictions  of  electro-optic 
properties  to  be  made  with  a  minimum  of  experimental  materials  characteri zation 
input . 

Finally,  the  experimental  aspect  of  this  program  has  been  quite  fruit¬ 
ful.  Original  characteri zation  measurements  of  the  linear  electro-optic  effect 
has  been  completed  in  TeC^,  and  begun  on  another  compound:  Tl^AsSe-^.  The 
studies  on  this  latter  compound  should  be  complete  by  late  1931. 

1.4  Tec  hnica  1 Issues 

The  most  significant  technical  issue  still  to  be  resolved  is  the  iden¬ 
tification  of  those  generalized  aspects  of  atomic  configuration  and  chemical 
bonds  that  will  allow  us  to  calculate  the  requisite  electronic  and  lattice 
responses  to  external  perturbations  (electric  fields  and  strains). 

1.5  Report  Summary 

The  remainder  of  this  report  contains  the  details  of  our  technical 
approach  and  the  results  for  both  the  theoretical  and  experimental  aspects  of 
this  program.  Finally,  there  is  a  short  summary  of  progress  and  some  of  our 
future  efforts.  Appendix  A  is  the  paper  on  the  experimental  efforts  on  TeO^* 
submitted  to  the  Journal  of  Applied  Physics. 
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2.0.  TECHNICAL  APPROACH 


2.1  Theoretical 


2.1.1  Electro-optic  Coefficients 

A  requirement  of  any  theoretical  endeavor  is  that  the  theory  must 
include  a  connection  to  the  experiment.  In  the  case  of  the  electro-optic 
effect,  the  measurement  gives  an  electro-optic  coefficient,  r^j,  which  is 
defined  by  the  relation^'® 


rijkEk 


or  in  reduced  notation 


(1) 


AB. 


(2) 


where  aB^j  is  the  change  in  the  relative  optical  dielectric  impermeability  and 
Ej<  is  the  applied  electric  field.  When  the  frequency  of  the  perturbing  electric 
field  is  greater  than  the  piezoelectric  response,  the  crystal  is  considered  to 
be  "clamped"  and  the  measurement  gives  a  constant  strain  electro-optic 

C 

coefficient,  r....  Elastic  deformation  occurs  for  a  low  frequency  (or  dc) 

1  j  K 

perturbing  field  and  this  "undamped"  situation  leads  to  a  constant  stress 
electro-optic  coefficient,  rj...  The  impermeabil ity  is  defined  by 

1  J  K 


B 


8E 

ij  "  M) 


f  =  )ij  '  (nr)ij 


ij 


(3) 


4 

C3485A/jbs 


Rockwell  International 

Science  Center 


SC5266.1FK 

where  is  a  dielectric  constant,  n  refers  to  a  refractive  index  and  is 
the  inverse  dielectric  tensor. 

The  change  in  impenneabi  ]  i  ty  with  applied  field  can  be  easily  relate’, 
to  a  change  in  refractive  index.  Let  n^  be  the  effective  refractive  index  fjr  a 
particular  set  of  propagation  and  polarization  directions  with  no  applied 
electric  field.  Then  let  n ^  he  the  effective  index  for  the  same  set  of 
propagation  and  polarization  directions  but  now  with  the  electric  field 
applied.  The  change  in  impermeability  is 


S3  = 


nj  -  "'2  (nj  -  h^){n1  +  n2) 


n< 


nfn2 


n  1  n*2 


( 4  ’ 


and,  in  the  approximation  that  -  n  and  .n^  -  =  -An,  it  becomes 


aB  =  -2tn/n3 


(5) 


The  same  result  can  be  achieved  from  differentiation: 


dB  d  ,1  , 
3n  ‘  '•  n* 1 


(6) 


The  explicit  relationships  between  An  and  aB^  have  been  derived^  for  optically 
isotropic,  uniaxial  and  biaxial  crystals  with  the  light  propagating  in  an 
arbitrary  direction.  For  example,  in  uniaxial  crystals,  the  change  in  the 
ordinary  refractive  index  is 


-n. 


Ano  = 


So  aB..  +  Sf  aB 


'll 


I* 


1  au22 


2S1S2aB12 


TT 
* 2 


(7) 
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for  a  propagation  direction  given  by  unit  vector  $  =  *  (0,0.1).  The 

expression  for  the  change  in  extraordi nary  refractive  can  be  written 


.  — T~- 


-2c?. 


c2e2. 


SjS,a9jj  + 


f  S1*S?  ‘  ?a°3'5 


?YzY 

i  ? 


aR 


1? 


-  2SzVYV 


aR 


?R 


- 


1? 


where  n  (eff)  is  the  effective  extraordinary  refractive  index  for  the 
propagation  direction  $  =  '  ,S?  ,S^  : 


n0(eff) 


n  n 
o  e 


/n7^T'ST 


Relating  the  impermeability,  and  hence  the  electro-optic  coefficient,  to  a 
change  in  refractive  index  is  of  fundamental  importance  because  1)  experi¬ 
mentally,  onc  nominally  meas  ires  this  change  in  the  index  of  refraction  as  a 
function  oe  applied  electric  field  and  2)  theoretically,  one  calculates  the 
perturbation  of  the  refractive  index  due  to  the  applied  field. 

2.1.2  Suscepti bi 1 i ty 

Traditional ly,  the  refractive  index  is  related  to  the  microscopic 
nature  of  a  crystalline  solid  via  the  polarizability  and  susceptibility.^  In  CGS 
units 


n  =  /I  +  4n  x  (10) 

for  isotropic  materials  in  the  linear  regime.  The  electro-optic  effect  can  be 
viewed  as  one  of  many  nonlinear  perturbations  to  this  linear  susrepti bi 1 ity;  in 
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particular,  it  corresponds  to  a  second-order  susceptibility.  There  are  many 
conventions  used  for  representing  the  nonlinear  susceptibility  coefficients. 
One  method,  which  is  instructive  but  not  rigorous,  is  to  consider  a  Taylor's 
series  expansion  of  the  induced  pol ari zat ion  ,  P,  in  terns  of  the  three  compo¬ 
nents  of  the  total  electric  field,  , 


o 

i 


4  aij 
J  J 


E'  +  4  )  b.  + 

J  2  K  1JK  j  x 


(ID 


where  the  expansion  has  been  done  about  the  point  t"  ~  0  and  it  has  been  assumed 
that  the  induced  polarization  for  zero  electric  field  is  zero,  i.e.,  P..  (I:  =  0) 

=  0.  Note  that  the  expansion  coefficients  are  a^,  b^,  etc.  Now  let  the 
total  electric  field  be  a  superposition  of  two  electric  fields  at  frequencies  ^ 
and  Lo  ^ 


*t 

t 


r;  +  f  •, 


(12) 


In  component  form. 


=  E1*’1  +  E^'- 
m  m  m 


Substituting  Eq.  (13)  into  Eq.  (11)  yields 


(13) 
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■jau!Ej‘*E“i,t?|kbiWEj“Ej2)(Er‘Ei:2!* 
•  jaijrjr‘  *  f  uEf  *■  jk-i-  *  jk-f-  E?Er 
*  i. -’f  *  4,^  Ek,Ej2  *  - 

J  K  J  K 


(14) 


Since,  in  the  summations  over  j  and  k,  j  and  k  are  only  "dummy"  indices,  it  is 
possible  to  interchange  j  and  k  in  the  last  temn  and  then  conbine  the  last  two 
terms  by  factoring  out  the  electric  field 


b  .  .  b  .  . 

P  =  a  FIJi  +  ^  a  F  J2  +  ■  _Ll<.  r  J l F  J a  +  ■■  — LI.  f^F"*2 

’  J  'J  i  j  m  2  Ei  Ek  Ik  2  EJ  Ek 


*ikKk*b,kj'ErEr  • 


(15) 


Now,  explicitly  include  the  frequency  dependence  of  the  fields: 


Ewi  =  cos  co,  t 
m  m  1 

and 

E“2  =  E^  cos  u>,t  . 
mm  2 

Substituting  Eqs.  (16)  and  (17)  into  Eq.  (15)  gives 


(16) 


(17) 
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b, 


Pi 


COS  uijt 


I  Vj1’  “S  “l1  *  \  Vf’  “s  V  *  E'1)Ek') 

J  J  J  K 

i  -4s-  Ej2M2)  co$2  "2*  * ),  i  *  W  E'”Ek2)  cos  "i1  cos  v 


J* 


Jk 


Using  the  trigonometric  relations  of 

cos  a  cos  ?,  =  -  cos  +  s)  +  cos  {.x  -  ii)  a” b  j  cos  (2  0  +  1  , 

P1  *  \  a1j  Ej l)  cos  "l1  +  l  a1jEj2)  C0S  “2*  +  l  ~^T  E^M^fcos  2<„  ,t  ♦  1) 

J  j  jk  J  1 

+  Er’E(k)(C0S  V  +  +  lk  7rbijk +  «»  (-nj  ♦  .2)t  ,  13) 

*  ]k  T  fhijk  +  bikj^jI>Ek?)  cos  (-!  -  «2)t  . 

Expressing  Pi  as  a  superposition  of  its  various  frequency  components,  i.e. 


p.  =  p°  +  p';’i  +  p'.J2  +  p2\  +  p2j.-  +  p'.Ji_J2  +  p:’i+  v 

iiiii  ii  i  ’ 


(19) 


these  separate  components  can  be  written  out  as 


p°  = 

'  JJJ.  f E ( ^ E f ^  +  E^E 

1 

jk  4  ^  k  J 

k  1 

p^i  = 

Mi 

)  a  .  .  E  i  '  cos  uj,t  , 

pw2  = 

' \A2) 

cos 

U)o  t  § 

i 

j  ij  J  1 

i 

3  n  J 

2 

p2“i  = 

Mp-  Ej1)£kl)  “5  V 

9 

pfJ2  ■  l 

Jk 

bijk 

4' 

E^E^  cos  2ui2t 

pWj-U)2  _ 

lkKjk  *  iWe(1,e(2) 

COS 

(li)j  - 

» 

and 

pW[  +u>2  _ 

lkKjk  *  W^'M2’ 

cos 

• 

(20) 
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By  examining  the  frequency  dependences  of  each  of  these  terms,  there  is  an 
obvious  correspondence  to  the  physical  phenomena  of  optical  rectification, 
second  harmonic  generation,  and  sum/difference  mixing.  In  addition,  for  the 
special  case  of  ^  =  0,  the  sum  and  difference  mixing  terms  are  equal  and 
correspond  to  the  linear  electro-optic  (Pockel's)  effect,  which  is  an  induced 
polarization  at  the  same  frequency  as  the  optical  field  and  is  proportional  to 
the  dc  electric  field. 

To  define  the  higher  order  susceptibilities,  one  can  write  the  induced 
polarization  in  terms  of  increasing  powers  of  electric  field:^ 


P; 


E. 

J 


} 

h 


E.E. 
J  < 


jks 


(3) 

(i  jkt 


EjV;c 


(?D 


where  x(*)  Is  the  linear  or  first-order  susceptibility,  x(^)  is  the  second-order 
susceptibility,  etc.  Within  this  context,  the  linear  Taylor  series  expansion 
coefficients  are  identical  to  the  first-order  susceptibility  tensor  elements. 
Combinations  of  the  second-order  Taylor  series  expansion  coefficients  correspond 
to  the  second-order  susceptibilities.  While  the  above  procedure  is  not 
completely  rigorous,  it  does  provide  the  appropriate  correspondences. 

An  alternate,  but  somewhat  similar,  approach  in  specifying  suscepti¬ 
bilities  is  to  define  the  ith  component  of  the  polarization  density  vector  as 
the  negative  of  the  partial  derivative  of  the  total  energy  density  with  respect 
to  the  ith  component  of  electric  field.  Then,  by  expressina  the  total  energy 
density  as  a  Taylor  series  expansion  in  the  electric  field,  the  susceptibility 
coefficients  can  be  related  to  the  derivatives  of  this  total  energy  density  with 
respect  to  the  various  fields.  Still  another  method  has  been  invoked  in  defin- 
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ing  the  second-order  susceptibility,  which  explicitly  emphasizes  the  frequency 
dependence.  A  general  form  of  the  second-order  polarization  can  be  written^ 


P ,  (r^)  -  ;  ;  (-■>.  ui* .  W)  £.(?,,•)  Ek(flU-u')  d^' 


jk 


(22) 


■+  ★  y 

where  negative  frequencies  are  permitted  by  assuming  E]  (r,J)  =  E.(r,-u>). 

The  specific  definitions  for  three  second-order  susceptibility  are :  ^ 


i) 

Second  Harmonic  Generation 

p2j 

*  '  v?:’  C 

jk  *1Jk  J  K 

( 23 ; 

ii) 

dc  Effect  (Optical  Rectification) 

P°  * 

0  r  jO  p  ’  J 

x,ji  j  k 

(24) 

iii) 

Electro-Optic  Effect 

py- 

V  U)  rO  r'j) 

jk  njk  Ej  Ek 

(25) 

Note  that  these  second-order  susceptibilities  ane  assumed  to  all  be  distinct 
tensor  quantities.  The  symmetry  in  the  subscript  indices  for  the  three  types  of 
susceptibilities  is  as  follows.  Second  harmonic  generation  and  optical  rectifi¬ 
cation  susceptibilities  are  symmetric  in  j  and  k  since  interchanging  two  optical 
fields  of  the  same  frequency  is  physically  insignificant.  The  electro-optic 
susceptibility  is  symmetric  in  i  and  k  as  seen  from  a  quantum  mechanical  formu¬ 
lation,  utilizing  the  electric  dipole  approximation,  or  as  seen  from  a  thermo¬ 
dynamics  argument.  This  latter  argument  can  be  understood  by  investigating  the 
following  question:  "What  is  the  total  susceptibility  that  the  optical  field 
'sees'  when  the  dc  field  perturbs  the  first  order  susceptibility?"  The 
polarization  is 
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pi  ■  i  us;1  *  i  x?  °k  e  j )  t 


(26) 


where  x ^ j k  in  Ecl*  (25).  For  a  principal  axes  coordinate  system, 

x\lJ  is  diagonal.  To  simplify  the  total  susceptibility,  consider  the  dc  field 
to  be  along  x  i.e.,  E°  =  E°  and  E°  -  E^  =  0.  Then 


lel\ 

/[  (!)  +  eo  Eo1 

[  Lm  xm  xJ 

eo  c  o 
x  112E  x 

eo  .  o 
x  113E  x 

M 

= 

eo  ro 

X21 1  Ex 

r  (1)  eo  0-, 

LX22  x212  xJ 

eo  ro 
x213  Ex 

E2 

Vi 

I  eo  o 

(  X3I i  Ex 

eo  to 

X312  bx 

r  Cl)  .  eo 
Lx33  X31 3 

0, 

\^l 

(27) 


In  order  that  this  total  susceptibility  tensor  be  symmetric,  which  is  a 
requirement  of  thermodynamics ,  it  is  necessary  that: 


eo  _  eo 
xijk  xkji 


(28) 


i.e.,  the  "electro-optic"  susceptibility  is  symmetric  in  the  first  and  third 
indices.  This  implies  that  the  "condensed  notation"  is  obtained  by  contracting 
the  two  indices  which  do  not  refer  to  the  dc  field.  In  contrast,  for  both 
second  harmonic  generation  and  optical  rectification,  the  two  indices  of  the 
optical  fields  are  "condensed." 

Finally,  by  equating  the  change  in  refractive  index,  due  to  the 
electro-optic  susceptibility  perturbing  the  linear  susceptibility,  with  the 
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change  in  impermeability  via  the  electro-optic  coefficient  (see  Eq.  (5)),  a 
direct  relation  between  the  second-order,  electro-optic  susceptibility, 

and  the  electro-optic  coefficient,  r ^ j k ,  may  be  obtained.  A  simple 
technique  for  approximating  this  relation  is  to  write  the  perturbed  refractive 
index  as 


n  +  on  = 


V‘ 


*  ♦  „me! 


1  +  4n 


;(1)  i 

'  1  +  4irxU', 


i  1/2 


\  4nxe°E 
=  n  1  +  -Jkr- 


;?9) 


and  then  use  the  binomial  expansion  on  the  square  root  to  give 
n  +  An  =  n  '1  +  ^-  [ l~r~nT~ 3  +  •  •  •  )  • 


Then,  by  referring  to  Eq.  (5), 


An  -  —  -  -j— 


which  implies  that 


-  -  -P- 


and,  since  aB  -  rE  from  Eq.  (1),  one  obtains 


r  -4  tty60 

r  -  -4-  . 


(30) 


(31) 


(32) 


(33) 


Performing  the  rigorous  matrix  inversion  required  to  convert  the  full 
susceptibility  tensor  to  the  impermeability  tensor  yields  the  relation7, 
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-4ux 


r. 


ikj 


eo 
cii bkk 


(34) 


where  =  n?  =  1  +  4nx-^  the  optical  dielectric  constant.  Note  that  con¬ 
vention  implies  "condensing"  the  first  and  second  indices  in  r^j  while  "condens- 
ing"  the  first  and  third  indices  in  x.  ...  However,  in  both  cases,  the  "con- 
densed"  indices  do  not  refer  to  the  static  electric  field.  Also,  the  dielectric 
constants  reflect  directions  other  than  that  specified  by  the  static  field. 

Let's  verify  Eq.  (34)  for  a  specific  example.  Consider  symmetry  point 
group  422  with  the  dc  electric  field  applied  in  the  <100>  direction.  Then,  from 
combining  Eqs.  (2)  and  (3)  and  using  the  electro-optic  matrices  in  Kaminow  6 
Turner,7  the  perturbed  impermeabi 1 i ty  is 


B  = 


h 

0 

0 1 

1/n0 

h 

0 

0 

l/n| 

0 

0 

0 

0 

+ 

r4i 

0 

0 

0  l 

0 

-r41 

0 

i 0 

0 

°/ 

/E*\ 

0 

\o  / 


/  lK 
0 

0 


1  /  n  0 
r41Ex 


r41Ex 

l/ne  / 


(35) 


where  nQ  is  the  ordinary  index  and  n0  is  the  extraordinary  index.  Referring  to 
Eq.  (26),  the  perturbed  susceptibility  tensor  can  be  written  as 


x  = 


/  Xo 
0 
0 

\ 


0 

Xo 

eo 


x312Ex 


eo  r 

x213Ex 
Xe 


(36) 
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since  x^  =  X31?  -  =  ‘x321  are  the  on^  nonzero  second-order  susceptibi¬ 
lities.  Note  the  Xo  ancl  Xe  dre  the  iinear  ordinary  and  extraordinary  suscepti¬ 
bilities,  i.e,  n2  =  1  +  and  n2  =  1  +  4-nXe.  From  Eqs.  (3)  and  (10), 


1  0  0 


0  1  0  +  4n 


0  0  1 


0  0 


u  x0  *X 
0  xe 


n?  0  0 

0 

0  n2  4ttAx 
0  4hAx  ne 


GO  GO 

where  ax  bas  been  defined  as  ax  -  X2i3^x  =  X312'x*  ^sin9  tlle  Gauss-Jordar, 
matrix  inversion  technique  to  invert  8"^  gives 


1  .  (4kAv)? 

nT  n^-r-fexT 


-4ttAx 


-4tiAx 


-4itAv 
n  2 
0  e 


0  44*  K 


in  the  limit  that  ax  is  small  so  that  terms  containing  (AX)n  where  n  >  1  are 
taken  to  be  zero.  Now,  equating  the  matrix  elements  in  Eqs.  (35)  and  (38) 
implies  that 


r  F  =  -41 & 

r4i  x  n£n Y 
0  e 
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_4tIx213  _  '4tix312 

r41  "  ~n7nT  'rv^iP 

o  e  o  e 

which  is  in  exact  agreement  with  Eq.  (34). 


(39) 


2.1.3  Units 


A  recurring  problen  encountered  when  dealing  with  various  nonlinear 
quantities  is  comparing  and  converting  between  MKS  and  CGS  units.  It  therefore 
seems  appropriate  to  summarize  some  of  the  conventions  used  in  presenting 
experimental  data,  along  with  a  few  conversion  factors. 


K 


The  experimental  electro-optic  coefficient,  r is  generally  expressed 
in  the  MKS  units  of  (meter/volt),  whereas  the  second  harmonic  generation  coeffi¬ 
cient,  di j  ,  has  often  been  given  in  the  "catch-all"  units  of  (esu)  and  in  the 
MKS  units  of  (meter/volt)  or  (cool omb/vol t^) . 

The  ambiguity  of  tne  "esu"  can  be  cleared  up  by  explicity  expressing 
the  esu  in  CGS  units.  This  ambiguity  stems  from  the  multiple  definitions  for  an 
"esu."  That  is,  1  esu  =  1  statcoulomb  for  charge,  1  esu  =  1  statvolt  for  poten¬ 
tial,  1  esu  =  1  statamp  for  current,  etc.  For  second  harmonic  generation 
coefficients  (or  any  other  second-order  nonlinear  coefficient),  an  "esu"  in  CGS 
units  is  defined  as 
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=  8.854  x  10 


-12  _ coulomb  _ 

newton-meter7 


.  newton-meter 
x  [ - voTt  = 


8.854  x  10 


-12  coulomb 
Volt -meter 


(41) 


for  several  different  reasons.  First,  it  turns  out  that  e0  is  dimensionless  and 
numerically  equal  to  (4 tt ) ~ ^  when  expressed  in  CGS  units: 


8.84  x 


in~ 12  coulomb 
'  vo ft -meter 


2.998  x  104  statcoulomb 

"  2-998  .1 0 : l .vol  ts' 

coulomb 

statvoft 

X 


_  1  meter 

100* ’cent  Cne  ter 


#- 

1_ 

- 


statvolt  -  centimeter 
statcoufomb 


(42) 


Secondly,  there  are  two  conventions  for  the  relation  defining  the  second  harmonic 
generation  coefficient.  These  two  relations,  referring  back  to  Eq.  (23)  for 
comparison,  are 


and 


E 


,o 

j 


E 


j 

k 


}kc0dijk  Ej  Ek 


(43) 

(44) 


where  the  two  indicated  second  harmonic  generation  coefficients  differ  by  a  factor 
of  eQ.  That  is  to  say, 


, coulomb » 

(ToTt^'J 


Go  dijk 


rmeter . 

f~oTfJ 


(45) 


Note  that  other  second-order  nonlinear  coefficients,  such  as  the  electro-optic 
coefficient,  which  can  be  expressed  in  units  of  (meters/volt),  usually  do  not 
require  this  factor  of  e0  when  converting  to  other  units.  (Instead,  use  the 
conversion  factor  (4ne0)  which  is  dimensionless  and  equal  to  unity.) 
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As  a  specific  numerical  example,  consider  the  results  of  two  different 
second  harmonic  generation  experiments  on  TeC^.  The  first  experiment^  used 
Eq.  (44)  to  define  d^  and  reported  a  value  of  |d^|  =  0.69  *  10"^  m/V.  The 
second  experiment^  used  Eq.  (43)  to  define  and  the  result  was  |di^|  =  1.45  x 
10“^  esu.  In  order  to  directly  compare  these  two  values,  use  Eq.  (45)  on  0.69  x 
10'1"  m/V  value  so  that  the  same  definitions  for  the  coefficient,  d^,  is  employed 
and  then  convert  to  esu: 


14' 


=  e0  (0.69 


10'12  m/V)  =  6.1  X 


10 


-24  coulomb 
voftf2” 


2^993  x  10  statcoulomb 
'(fo'Jfomb  ~ 


2.993  _x_l 0 2  volt 

2 

1  esu 

TfatvoTt” 

s  t  a  t  c  o  uTomb /s  t  a  t  v  o  1 1 

(46) 


which  is  in  reasonable  agreement  with  the  other  value  of  1.45  x  10_C)  esu. 

The  electro-optic  coefficient,  r- ,  the  second  harmonic  generation 

(2) 

coefficient,  d^j,  and  the  second-order  susceptibility,  X > . ' ,  all  have  the  same 
dimensions,  as  is  evident  from  Eqs.  (25),  (34),  (43),  and  (44).  Provided  that 
these  coefficients  are  consistently  defined,  in  contrast  to  the  situation  in 
Eqs.  (43)  -  (46),  a  particular  coefficient  may  be  expressed  in  units  of 
(meter/volt),  (coulomb/volt2),  (esu),  etc.  by  using  the  following  conversions: 

1  meter  A  0  ocn  iri-12 

"volt —  *  4tt  x  o«oo4  x  10 

1  coulomb  [2.998  x  10q  statcoulomb  |  |  2.998  x  102  volt~|  f _ 1  esu _ 

voTt2  x  L  coulomb"  “  j  x  L  statvoTt”  J  x  [_ sta t cVuTomb/stat voTt7 

=  2.695  x  10+14  esu  (48) 


coul omb  .  . .  , 
volt -meter J  ~  i,11J 


10 


-10 


coul omb 
“To!  t7  “ 


(47) 
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2. 998  x  10“  esu 
(49) 


1  meter 

100  centimeter 

2.998  x  10'  volt 

vol  t 

meter 

X 

'statvdTt  . 

iriu  centimeter 
statvdTt ~ 

(50) 


where  each  quantity  in  the  square  brackets,  [  ],  is  dimensionless  and  has  a 
value  of  unity.  As  a  numerical  example,  utilizing  these  conversion  factors,  let 
us  determine  the  electro-optic  susceptibil ity  (in  esu)  for  GaAs  associated  with 
an  experimental  electro-optic  coefficient  of  r^  =  1.2  x  10‘12  m/V  and  a 
refractive  index  of  n  =  3.60.  Using  Eqs.  (34)  and  (49), 


•'  2h  ’  -(3.60)4  (1.2  x  *1J'12  m/V )  2. 993  x  104  =  -4.8  x  10'7  esu. 

(51) 

Finally,  f rom  a  theoretical  point  of  view,  the  most  convenient  repre¬ 
sentation  for  these  dimensions  (see  ROM  section  below)  is  charge^/energy^.  To 
be  compatible  with  this  form  (and,  in  addition,  combining  the  results  from 
Eqs.  (49)  and  (50)),  Eq.  (40)  can  be  rewritten  as: 


,  _  1  statcoul  _  1  centimeter  _  1  statcoulomoj 

u  ~  statVoTt 7~  ”  statvoff  erg-  " 


(52) 


since  1  statcoulomb  =  1  centimeter-statvolt  and  1  statvolt  =  1  erg/statcoulomb. 
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2 . 1 . 4 . 1  Introduction 

One  form  of  tight-binding  theory,  the  Bond  Orbital  Model  (BOM),  has 
previously  been  used  to  describe  numerous  electronic  and  lattice  properties  o' 
crystalline  solids.  As  originally  formulated  by  Prof.  W.A.  Harrison  of  Stanford 
University,  the  BOM  was  applied  to  tetrahedral ly  coordinated  solids,  in  parti- 
cular,  to  the  elemental  and  zincblende  semiconductors.  *-  Specifically,  this 
model  has  been  used  to  calculate  the  dielectric  properties  of  tetrahedral 
solids,  including  the  refractive  index  for  the  zincblendes*  and  for  the 
chal copyri tes  along  with  second-order  optical  susceptibilities  for  wortzite 
and  zincblende  compounds.^7  Further  refinements  in  the  tight-binding  parameter¬ 
ization  have  previously  occurred^’19  and  are  still  in  progress.^0  These  new 
sets  of  parameters  do  not  affect  the  basic  tenets  of  the  BOM  formalism  itself, 
but  enter  only  when  applying  the  theory  to  certain  specific  cases.  The  BOM  has 
recently  been  reviewed  in  detail  and  its  scope  of  applicability  has  been 
extended  to  include  other  crystal  symmetries. 

2 . 1 . 4 . 2  Susceptibil ity 

The  lowest-order  susceptibility  can  be  expressed  with  the  BOM  through 
an  expansion  of  the  total  energy  in  the  electric  field.  The  second-order-  sus¬ 
ceptibility  is  then  realized  by  carrying  the  expansion  to  one  higher  order  in 
the  electric  field.  The  simplifying  concept  of  the  BOM,  which  permits  a  "first 
principles"  prediction  of  the  dielectric  properties,  is  the  following.  By 
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calculating  the  change  in  polarization  in  each  "chemical"  bond  caused  by  the 
optical  electrical  field,  one  obtains  the  total  susceptibility  by  summing  up  th- 
contributions  from  each  individual  bond.  When  the  crystal  geometry  (i.e.,  each 
bond)  is  perturbed  by  a  dc  electric  field,  the  polarization  associated  with  each 
bond  changes.  Consequently,  the  susceptibility  is  also  modified  and  it  is  this 
change  in  susceptibility  with  respect  to  the  dc  electric  field  which  corresponds 
to  the  linear  electro-optic  (Pockels)  effect. 

An  alternative,  but  equivalent,  procedure  is  to  directly  associate  the 
change  in  polarization  per  bond  due  to  both  the  optical  and  dc  electric  fields 
wi to  the  second-order  susceptibility.  In  the  BOM,  the  dipole  moment ,  p,  of  each 
bool  is  written  as 

P  *  -eyiap  =  -  9Y3  V3(V?  m  ^)*1/2  (53) 

where  e  is  the  electronic  charge,  r  is  the  "center-of-gravi ty"  for  hybrid  orbi¬ 
tals,  3  is  the  interatomic  distance  vector  between  the  two  bonding  atoms,  a,,  is 
the  "polarity"  of  the  bond,  V?  is  the  "covalent"  bond  energy,  and  V3  is  the 
"polar"  bond  energy.  With  reference  to  Eq.  (22),  the  first-order  susceptibil¬ 
ity*  a  ^  the  average  change  in  dipole  moment  responding  to  the  optical 

electric  field,  £0(3t.  In  addition,  the  second-order  susceptibility,  »  is 
the  average  change  in  dipole  moment  responding  to  both  and  the  dc  electric 
field,  £dc.  These  two  susceptibilities  can  be  expressed  as  derivatives  in  these 
electric  fields  by 


,i'i  .i  E  J-L,- 

v  bf^s  dEopt 
and 


(54) 
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(55) 

where  the  average  has  been  obtained  by  systematical  ly  summing  over  each  type  nr 
bond  in  the  crystal  and  then  dividing  by  the  volume,  v,  occupied  by  those  bonds. 

The  derivative  of  the  dipole  moment  per  bond  with  respect  to  one  or 
both  of  the  electric  fields  can  be  written  in  terns  of  partial  derivatives  of 
the  dipole  moment  per  bond  with  respect  to  various  BOM  parameters  or  interatmic 
distance,  multiplied  by  the  change  in  that  particular  BOM  parameter  as  a  func¬ 
tion  of  the  appropriate  electric  field.  This  is  just  an  application  of  the 
chain-rule  for  partial  di f ferenti ation. 

For  example,  the  derivative  in  the  first-order  susceptibility  becomes 


(2) 


.1  y  _  j{p 


v  bonds  dE0ptdEdc 


f  dl__'  -  '  ld  +  (  1V2  n  +  '  w3_p_^ 

dE0pt  -  L  3F.°Pt  ' ad  Vpt),'9V 


(55) 


Each  partial  derivative  can  be  evaluated  separately.  By  making  the  assumption 
that  the  frequency  of  the  optical  electric  field  is  high  enough  so  that  the 
lattice  cannot  respond  to  it  (i.e.,  greater  than  the  Reststrahl  frequency),  the 
interatomic  distance  can  be  taken  to  be  independent  of  the  optical  electric 
field.  In  other  words, 


3d 


r*-l  =  0 


3E 


opt 


(57) 


Since  the  covalent  energy  can  be  represented  as  being  proportional  to  the  square 
of  the  interatomic  distance,  i.e., 
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2.16  *■ 
~cP"~ 


(58) 


where  1i  is  Planck's  constant  divided  by  2n  and  m  is  the  electron  mass  (li^/m  = 
7.62  eV-£-)»  the  chain-rule  gives 


Therefore,  only  the  "polar"  bond  energy  term  in  Eg.  (56)  is  nonzero.  Evaluating 
this  tern  yields 


and 


=  -  ey(t  V'  ■  Vv) 


•  E0pt./2 


(60) 


(M) 


where  E0,,t  is  the  unit  vector  giving  the  direction  of  the  optical  electric 
field.  This  last  equation  results  from  the  change  in  energy  due  to  the  inter¬ 
action  of  the  dipole  with  the  total  electric  field.  This  dipole  energy  is 
-  ye3  •  (£0pt  +  fdC)/2.  Finally,  by  substituting  Eqs.  (57),  (59),  (60)  and  (61) 
into  Eq.  (56),  the  first  order  susceptibility  can  now  be  written  as 


in  .  *  T 

^  ^  bonds  (V|  + 


(62) 


where  E°pt  is  a  unit  vector  in  the  direction  which  represents  the  optical 
electric  field. 
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Followimj  a  similar  procedure  for  the  second-order  susceptibility,  the 
derivative  of  Eq.  (56)  with  respect  to  the  dc  electric  field  becomes 


f-4SUl  •  ,r  >2-L-t  . 


u  p  \  -  f  J  i  If  i,'  o  ji  .  r  £■  w  o-p _ 

"  v37“pTTjF  jv,sa  Vc  rVv2 


where  any  terms  involving  the  variation  of  interatomic  distance  with  the  opti< 
electric  field  have  been  dropped.  Noting  that 


•1  »  -  ye3  .  t  /2 


,.»f  .  3e^V3_  ,« 

la'vI  ‘  (Tf.'Vp2  ’ 

in  addition  to  using  Eq.  (61),  implies  that  the  "electronic"  portion  of  the 
second-order  susceptibility  can  be  written  as 


xjjk  (electronic) 


Y^v|v3di;3.  eJc)(J  •  E°?t) 

*-  ^  - - - - V  TFT - ‘ - 


bonds 


(»i  *  V1 


This  contribution  to  the  second-order  susceptibility  has  been  called  the  "elec¬ 
tronic"  part  since  it  has  been  derived  from  the  deformation  of  the  electron  dis¬ 
tribution  in  the  bond.  The  deformation  of  the  electron  cloud  due  to  both  the 
optical  and  dc  electric  fields  results  in  a  change  in  the  "polarity"  of  the  bond, 

Another  contribution  to  the  second-order  susceptibility  arises  from 
those  terms  in  Eq.  (63)  which  involve  the  lattice  response  to  the  dc  field. 

This  "lattice"  contribution  has  been  further  divided  into  an  "ionic  displacive" 


% 
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part  and  "piezoelectric/photoelastic"  part.  The  "ionic  displacive"  component 
involves  only  the  relative  motions  of  the  atoms  whereas  the  "piezoelectric/ 
photoel asti c"  component  allows  a  piezoel ectric  strain  or  bulk  lattice  distor¬ 
tion.  The  reason  for  this  last  division  is  to  permit  compatibility  with  dis¬ 
tinction  of  a  "clamped"  or  "unclamped"  measurement  of  the  electro-optic  coef¬ 
ficient.  A  "clamped"  measurement  corresponds  to  a  second-order  susceptibility 
which  includes  only  the  electronic  and  the  ionic  displacive  lattice  contribu¬ 
tions,  and  not  the  piezoelectric/photoel astic  portion. 

The  "lattice"  contribution  to  the  second-order  susceptibility  is  ob¬ 
tained  by  combining  Eq.  (55)  with  the  first  two  terms  of  Eq.  (63).  Evaluating 
each  factor  in  this  latter  equation,  for  which  an  explicit,  simplifying  expres¬ 
sion  can  be  obtained,  implies 


JLiP_  =  .  eyy2  V-i  +  Vf 3/2 
av3act 


^  *  v 


-5/2 


and,  just  as  in  Eq.  (59), 


3V 

aE 


2  ,  _  ,  3d 


3V, 


ad 


dc  -  i vac)' 'var  -  “dv J l  "d~ 

aE  aE  aE 


-2V, 


(67) 

(6-1) 

(69) 


Then,  by  substituting  the  above  3  equations  along  with  Eq.  (61),  the  "lattice" 
contribution  to  the  second-order  susceptibility  can  be  written  as 
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*ijk(lattice)  =  fv 


y2H 


bonds  f V|  + 


Dd .  d.  . 

i V1  *  v!K-^nr)  *  2(v|  -  2v|)  r  (-’4-)  (J 

3lj  3Ej 


fopt 

Lk 


(70) 


where  a  distinction  is  made  between  the  magnitude  of  the  interatomic  distance 
vector,  d,  and  the  ith  component  of  the  interatomic  distance  vector  d;  (i.e., 

3  =  y'd.x.).  This  contribution  to  the  second-order  susceptibility  is  called 

M  1  1 

the  "lattice"  part  since  it  has  been  derived  from  atomic  displacements  (both 
relative  atomic  motions  and  bulk  piezoelectric  distortions)  due  to  the  dc  elec¬ 
tric  field.  In  fact,  the  lattice  response  to  the  dc  field  enters  into  the  for¬ 
mal  arization  as  a  single  factor,  (33/3Edc),  which  is  the  change  in  the  inter- 

J 

atomic  distance  vector  with  respect  to  the  dc  electric  field. 

In  summary,  the  second-order  susceptibility  has  been  described  in  terms 
of  the  Bond  Orbital  Model.  General  formulae  (Eqs.  62,  66  and  70)  have  been  de¬ 
rived,  which  permit  direct  calculation  of  the  first-order  susceptibility  and  the 
electro-optic,  second-order  susceptibility  for  any  compound  with  a  known  crystal 
structure.  The  necessary  BOM  parameters  (V£,  V3  and  y)  are  readily  available 
and  the  only  missing  ingredient  is  the  factor  that  quantifies  the  lattice  re¬ 
sponse  to  the  dc  electric  field.  Some  examples  on  evaluating  this  factor  are 
given  below.  The  second-order  susceptibility  then  can  be  converted  to  an 
electro-optic  coefficient,  r^j,  by  using  Eq.  (34)  in  conjunction  with  Eqs.  (3), 
(10),  and  (62). 

A  few  numerical  applications  for  the  first-order  susceptibility  and  the 
electronic  contribution  to  the  second-order  susceptibility  and  electro-optic 
tensor  are  illustrated  in  Table  1.  Comparisons  with  experiment  will  be 
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meaningful  only  when  the  lattice  contribution  to  the  second-order  susceptibility 
is  included  in  the  calculation  and  therefore  will  be  postponed  until  later. 


2 . 1 . 4 . 3  Appl icati on  to  Tetrahe dr  ally  Coordinated  Compounds 

The  prescription  for  calculating  susceptibilities,  given  in  the 
preceding  section,  has  been  applied  to  the  binary  zincblende  compounds.  Since 
all  bonds  are  equivalent,  the  summations  over  bonds  which  appear  in  Eqs.  (62), 
(66)  and  (70)  can  be  simplified  to  give  analytic  expressions  that  do  not  contain 
any  summations.  By  factoring  the  BOM  parameters  out  of  the  summation  because 
they  are  independent  of  bond,  the  summations  can  be  evaluate!  for  the  fou1' 


interatomic  distance  vectoi 

rs  : 

(1,1,1)  d//T 

(71) 

a'2)  = 

(1,-1, *1) 

d//3 

(72) 

a(3)  „ 

(-1 ,1,-D 

d//3 

(73) 

ii 

(-1  ,-1,1) 

d//3 

(74) 

where  d  =  /3  a/4  and  a  is  the  lattice  constant.  The  summation  in  Eq.  (62) 
becomes 

2  d,df  •  h2/ 3  ,or  1  ■  J  (75) 

4  bonds  1  J  (  0  for  i  t  j 

and  the  summation  in  Eq.  (66)  becomes 


2  VA  =| 

bonds  J  I 


+4/3  d3/9  for  i  *  j  *  k 

0  for  i  =  j  or  i  =  k  or  j  =  k 


(76) 


—  -  ■"*  *• 


1 
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However,  the  choice  of  signs  for  the  four  interatomic  distance  vectors  in 
Eqs.  (71)-(74)  was  arbitrary.  The  crystal  could  have  been  represented  equally 
well  by  =  (-1.-1.-1)  d//3,  <J(2)  =  (-1 ,1  ,l)d//T,  =  (l,-l,l)d//3  and 

3^  =  (1 ,1 ,-l)d//3.  Equation  (75)  is  invariant  to  this  sign  change  but 
Eq.  (76)  is  not;  the  right-hand-side  of  Eq.  (76)  must  be  negated.  Before  deal¬ 
ing  with  the  lattice  response  necessary  to  proceed  with  Eq.  (70),  the  simplifier 
expressions  for  Eqs.  (61)  and  (66)  will  now  be  written  as 


and 


c(1)  - 


(2) 


Me'  r2  V2  d2 
12(V2  +  V2)3"'77 


for  i  =  j 


for  i  *  j 


(77) 


xjjjj  (electronic)  =  -j 


/I  Ne3  V|V3  d3 

- 1— —  for  i  t  j  *  k 

24  (V2  +  V2)5^ 

0  for  i  =  j  or  i  =  k  or  j  =  k 


(78) 


respectively,  where  N  =  3/v  =  4/a3  =  S/"?/  16/d 3  is  the  electron  density  (i.e., 
the  number  of  electrons  in  four  bonds).  These  two  equations  agree  with  Eqs.  (4- 
28)  and  (5-13)  of  Reference  3. 

In  order  to  simplify  Eq.  (70),  the  latt’ce  response  in  the  binary  zinc- 
blendes,  induced  by  a  dc  electric  field,  must  first  be  estimated.  As  the  ini¬ 
tial  step,  consider  only  the  ionic  displacive  portion  of  the  lattice  response. 
Using  GaAs  as  a  specific  example,  the  relative  motion  of  a  Ga  cation,  in  a  dc 
field,  will  be  determined  with  respect  to  a  rigid  lattice  of  As  anions.  This 
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relative  motion  is  illustrated,  in  an  exagerated  fashion,  in  Fig.  1  for 
^dc  _  ^-dc  This  perturbed  interatomic  distance  vectors  can  be  written  as 


3(1)  *  ([1  -  5],  1,  l)d//3 

(79) 

3(2)  =  ([1  -  5]  -1,  -ljd/,'3 

(80) 

3(3)  =  '-[1  +  6],  1,  -1 )d//3 

(81) 

3(4)  =  (.[!  4  6],  -1,  ljd//l 

(82) 

where  u  =  ux  =  5d//3  x  is  the  relative  vector  displacement  of  the  Ga  cation.  By 
associating  two  lattice  vibrational  frequencies  with  a  "classical"  spring 
constant,  k,21 


K 


=  j 


80, 


) 

t 


where  C0  and  Cj  are  the  "bond-stretching"  and  "bond-bending"  force  constants, 
the  displacement  can  be  obtained.  The  force  of  the  displaced  Ga  cation  is 
equated  with  the  force  exerted  by  the  dc  electric  field  on  this  cation: 

t  =  Ku  =  ee*  (83) 

where  u  is  the  vector  displacement,  e  is  the  electronic  charge,  e^  is  the 
dimensionless  "transverse"  effective  charge  and  i:dc  is  the  applied  dc  electric 
field  corrected  for  the  dc  surface  charges.  Solving  for  the  displacement  gives 


3ee*  £dcd2 

5  =  TTc-r^T 


(84) 


and  e?,  C  and  C,  have  been  calculated  or  tabulated  in  Reference  3. 
T  o  1 
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Now,  referring  back  to  Eq.  (70)  along  with  Eqs.  (79)-(82)  and  (84), 
evaluate  the  lattice  response  factors  ( 3di /sEj c )  and  (ad/aE*?0)  for  each  of  the 
four  bonds.  The  first  factor  is 


^i_ 

3Ef 


■6d/(/J  eJc' 


E<ic 


-3ee^.  d2 

4TC"~8C'*y  for  1 


for  i  t  j 


*  J 


(85) 


which  leads  to  a  summation,  in  the  crystal  geometry  factors  in  Eq.  (70),  of 


£ 

4  bonds 


2  <k 
4  bonds 


for  i  =  j 


for  i  *  j 


=  +  4u?/EjC  for  i  =  j  --  k  (86) 

-  0 

V 

since  n  d.  -  -4u..  This  summation  is  taken  to  be  zero  when  retaining 

4  bonds  k  j  J 

terms  only  to  first-order  in  the  displacement.  To  first-order  in  5,  the  second 
factor  becomes 


ad  _  d/r±  2 5 /5  - 

~az  =  r«rc* - 

hj 


(87) 


where,  for  example,  the  -  sign  is  for  bonds  1  and  2  while  the  +  sign  is  for 
bonds  3  and  4  when  E^c  =  E  .  By  performing  a  binomial  expansion  on  the  squared 
root, 


32 

C3485A/jbs 


Rockwell  International 

Science  Center 


ad 


(t).  -4- 

1  '■>  -dc 


J 


3E( 


(±)j|u|  (±)j/J  ee^  d?- 

,'T  rd~c  4  rc0  + 


SC5266.1FK 


(88; 


where  |uj  is  the  magnitude  of  the  displacement.  Substituting  this  into  the  last 
term  of  Eq.  (70)  gives 


di  3d 

.  T“a  cT  “He 
4  bonds  aE. 

J 


fc\ 


>'3  ee^.  d 

=  Trrjn-c-)- 


E 

4  bonds 


(  +  )  -d-d, 

V  'j  1  K 


!  a- 


and  this  last  summation  reduces  to 


E 

4  bonds 


(+) . d. d. 

'  ;j  i  k 


j  -4d:/3  for  i  t  j  *  k 
I  0  for  i  =  j  or  i  =  t;  or  j  =  k 


(r . 


As  in  Eq.  (7b),  the  sign  is  dependent  upon  the  choice  of  sign  for  the  four 
interatomic  distance  vectors.  Finally,  by  replacing  the  results  of  Eqs.  (3b), 
(39),  and  (90)  into  Eq.  (70),  the  ionic  displacive  lattice  contribution  to  th-’ 
second-order  susceptibility  for  zincblendes  can  be  written  as 


xi?!  (lattice) 


/3  Ne3e*y2V|( -  2V| jd ^ 

24  IV2  +  W)5'2  (C  +  8C.) 
2  3  o  1 ' 


(91) 


where  N  is  again  the  electron  density.  Numerical  examples  using  this  last 
equation  along  with  Eqs.  (77)  and  (78)  are  presented  in  Table  2  for  the  simple 
tetrahedral  solids,  which  have  force  constants  listed  in  Table  9-1  and 
transverse  effective  charges  calculated  from  Eq.  (9-24)  of  Reference  3. 

The  BOM  susceptibility  formulae  for  the  zincblende  compounds,  i.e, 
Eqs.  (77),  (78)  and  (91),  can  be  readily  modified  to  accommodate  the  chalcopy- 
rite,  abc2»  crystal  structure.  In  the  "quasi -cubic"  approximation  the 


33 

C3485A/jbs 


Table  2 

Evaluation  of  the  Electro-Optic  Coefficients,  Including  the 
Ionic  Displacive  Contribution,  for  Some  Zinch lends 


Rockwell  International 

Science  Center 


SC  5266.  If- 


34 

C3485 A/ j  bs 


Rockwell  International 

Science  Center 


SCS266.1FP 


chal copyn te  is  regarded  as  being  a  zincblende  with  two  different  types  of  bonds. 
Then,  the  summations  over  bonds  in  Eqs.  (62),  (6b),  and  (70)  reduce  to  two  terns 
instead  of  a  single  term;  one  for  a-c  bonds  and  another  for  the  b-c  bonds.  Each 
of  these  two  terms  is  evaluated  separately  using  the  appropriate  set  of  'iOK 
parameters  and  interatomic  distance  which  characterize  the  corresponding  bond. 

For  example,  the  first-order  susceptibi 1 ity  of  Eq.  (77)  becomes 


(1) 


r2V?d2 

•rzV  2d2 

(V|  +  V2;3/2 

a-c 

b-c 

for  i 


0 


for  i  t  j 


The  second-order  susceptibilities  can  si  mi  1 a^ly  be  written  as  two  terns  with  the 
leading  numeric  factor  being  halved  to  compensate  for  the  separation  of  the 
contributions  from  the  two  types  of  bonds.  The  "quasi -cubic"  approximation 
leads  an  isotropic  first-order  susceptibility,  i.e.  x^  =  x^  =  Like¬ 
wise,  the  second-order  susceptibility  has  i  =  'x213  =  x312'  “ 

1  ?)  ( 2)  _ 

XH2  =  X ^ 3 1  *  ^  the  actual  symmetry  of  the  chalcopyrite  structure  (42m)  is 

employed  in  Eqs.  (62),  (66),  and  (70),  then  the  degeneracies  in  the  suscepti¬ 
bilities  are  broken:  x^j  =  x^  *  x^ 


r  (?)  .  (2), 

<132  "  *231  * 


<»  and  -,!?)  ■ .(?!) .  -,<?>  -  AJh  , 


*123  *321 1  '*213  *312' 


2 . 1 . 4 . 4  Lattice  Dynamics 

General  formulae,  which  can  be  used  to  determine  the  first-  and  second- 
order  susceptibility  tensors,  have  been  given  in  Eqs.  (62),  (66),  and  (70).  All 
three  of  these  formulae  require  a  knowledge  of  i)  the  crystal  geometry  which 
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specifies  the  interatomic  distance  vectors  and  ii)  the  BOM  parameters  that 
characterize  each  type  of  bond.  In  addition,  the  last  formulae  (Eq.  7'J) 
requires  a  knowledge  of  iii)  the  lattice  response  to  a  dc  electric  field.  Toe 
purpose  of  this  section  is  to  provide  one  detailed  procedure  for  evaluating  to  is 

,-j  r 

latte-  lattice  response  factor,  ’  pd/jE . 

[he  approach  taken  here  is  intended  to  give  only  the  ionic  displative 
contribution  to  the  lattice  response.  Hence,  the  calculated  electro-optic 
coefficients  should  be  compared  with  "clamped"  electro-optic  measurement. 

Briefly,  the  method  may  be  outlined  as  follows.  Consider  an  isolated 
a to n  (?Q),  which  is  free  to  move  in  response  to  a  "dc  1  electric  field, 
surrounded  by  a  rigid  lattice  of  N  nei ghbo-i ng  atoms  (4\,  2,  •••  N).  Let  atom 
#1  be  coupled  to  each  of  tuese  fixed  neighboring  atoms  by  a  classical  spring, 
which  is  characteri zed  by  the  sp-ing  constant,  ,  where  i  =  1,2,  •••  N.  For 
snaM  di spl acements ,  the  relative  motion  of  atom  #0  with  respect  to  the  rigid 
1  attic  can  be  determined  by  equatinj  the  force,  which  the  electric  field  exe-ts 
on  the  effective  charge  of  atom  #0,  with  the  force  that  the  N  springs  exert  on 
atom  #0  when  it  is  displaced  from  equilibrium.  The  spring  constants  are  quan¬ 
tified  by  treating  atom  #0  as  a  simple  harmonic  oscillator  and  solving  the 
equation  of  motion  for  frequencies  and  normal  modes  of  vibration.  By  associ¬ 
ating  the  frequencies  of  vibration  with  the  IR  absorption  (Reststrahl)  fre¬ 
quencies,  the  spring  constants  are  determined. 

In  order  to  specify  the  problem  in  more  detail,  let  f  be  the  vector  to 
the  equilibrium  position  of  atom  #0,  let  r  be  the  vector  to  the  displaced 
position  of  atom  #0,  let  u  be  the  vector  displacement,  and  let  i*.  be  the  vector 
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to  the  fixed  position  of  atom  #i  (i  =  1,  2,  •••  N),  as  illustrated  in  Fig.  2  for 
N  =  3.  Assuming  that  all  springs  are  unstretched  when  atom  #0  is  at  its 
equilibrium  position,  the  force,  ^  ,  on  atom  #0  due  to  the  ith  spring  is 


=  -k.' 

i  ^ 


ri 


r  - 


i  >  >  | 

r .  -  r 
1  1  o' 


>  -► 
|rV-’r| 


(93) 


where  j r .  -  r  |  is  the  unstretchei  length  of  the  i^  spring,  |r.  -  rj  is  the 
stretched  length,  ( j  r .  -  rj  -  |r^  -  r  |)  is  the  amount  of  stretch  in  the  i^*1 
spring,  and  (rv  -  r)/|r.  -  r|  is  the  direction  of  the  restoring  force  for  the 
ith  spring.  Note:  defining  the  restoring  force  as  being  directed  from  the 
displaced  position  of  atom  #0  to  the  fixed  position  of  atom  #1  assumes  that  the 
spring  "connections"  to  the  atoms  a^e  free  to  pivot. 


is 


When  the  dc  electric  field  is  zero,  the  equation  of  motion  for  atom  ai 


d2r 

3t7 


N 


ri 


n  -  if,  -  f0n 


-► 

ri 


(94) 


where  mQ  is  the  mass  of  a ton  #0.  The  simple  form  of  this  differential  equation 
results  from  the  assumption  of  an  isolated  atom  coupled  to  a  rigid  lattice. 

If  a  spring  which  connects  two  atoms,  both  "free"  to  respond  to  external 
influences,  is  included,  then  a  more  complicated  set  of  coupled  differential 
equations  must  be  solved.  The  displaced  position  of  atom  #0  is  related  to  the 
equilibrium  position  by 

•  ?  =  ?0  +  u  (95) 
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Fig.  2  Example  arrangement  of  "free”  atom  #0  coupled  to  three  "rigid"  atoms  by 
springs. 
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and  it  is  the  vector  displacement  u  which  undergoes  simple  harmonic  motion. 
That  is. 


u(t)  =  uQ  cos  (jt  +  $) 


and  therefore 


d2r 

d-:  r  0  d  <  'i  n 

dt?  =  dt7"  +  dt7  J 


u>~  J 


(96) 


(97) 


where  u  =  ax«  +  Ayy  +  a zz  is  a  small  quantity.  Then,  the  equation  of  motion 
become'-. 


vn  ,<?(r_p)+  K 

0  1  o'  1 


1  - 


lri 


!?i 


+  I 

n 


(r%  -  r)  =  0 


(9°. 


This  equation  represents  a  set  of  nonlinear  differential  equations  which  are 
coupled  in  the  three  unknowns,  ax,  Ay,  and  az. 

The  nonlinear  factor  in  Eq.  (98)  can  be  linearized  in  the  following  way. 
Let  <3io  =  -  ?Q  which  implies  that  dio  =  /Tx”~ * q )'z T Ty  ~ Y~V ~+  T2  [ ’z^T7. 
Then,  the  x-component  of  the  nonlinear  factor  in  Eq.  (98)  can  be  written  as 


i  -►  ■* 

r.  -  r 


— - - °—  (X.  - 

±1  1 


iri  -  ri 


di0(^i  -  X0  -  AX) 

x)  "  /T*""-  x0  -  AxP^Tyi  ”  y0  -  Ay77*+  {7.  -  zQ  -  wj* 

(xi  -  *0  -  Ax) 


1  -  (^T-)[(xi  -  xq)ax  +  (y.  -  yQ)Ay  +  (z.  -  Zq)az] 


10 


where  only  term  in  first  order  in  ax,  a y  and  az  have  been  retained.  By  using  the 
binomial  expansion  on  the  denominator,  this  nonlinear  factor  becomes 
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lfi  -  fol  , 

— - —  (x,  -  x) 

I  ±  ±1  1 


|r1  -  r | 


(X.  -  xQ  -  ax)  +  ^ 


TO 


(xi  -  xo)2ax  +  <xi  -  xo)(yi  -  yo)iy 
+  (xi  '  x0^fzi  '  Z0>AZ]  #  (9J' 


Finally,  the  linearized  equation  of  motion  can  be  written  as 


*  V2J  *  £  Iy-- -  f„)  •  4 


=  0 


(100) 


which  is  an  eigenvalue  equation  in  the  displacement  variable,  u. 


By  diagonl izing  the  matrix,  A,  whose  elements  are  gi  ven  by 


mn 


N  -K. 
=  D  — -1- 
i  =1 


-► 

ri  ‘ 


"72 


ri  - 


o  m 


r .  -  r 
l  o 


(101) 


the  eigenvalues,  m^2,  are  determined  and,  hence,  the  three  frequencies  are 
associated  with  the  three  normal  modes  of  vibration  for  a  given  set  of  spring 
constants,  (i  =  1,2,  •••  N).  An  appropriate  set  of  spring  constants  can  bo 
found  numerically  by  using  iteration  techniques,  in  which  the  K^'s  are  varied 
the  eigenfrequencies  obtained  match  the  experimentally  observed  lattice 
optical  phonon  frequencies. 


Once  the  spring  constants  are  known,  the  "dc"  electric  field  can  be  ap¬ 
plied  and  the  force  exerted  by  this  electric  field  will  exactly  balance  the 
forces  of  the  springs  on  atom  #0.  Specifically,  again  linearizing  the  forces  of 
the  springs  (see  Eq.  (99)),  we  obtain 


fdc  =  +  it  nrrri*  ^  -  **){{**  -  r+°}  * 


(102) 
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where  e^;  is  the  effective  transverse  charge.  Solving  this  equation  for  the 
di spl acement  vector,  u,  will  result  in  the  determination  of  the  ionic  displacive 
lattice  response  factor,  ( C ) . 

Initial  steps  have  been  taken  in  applying  the  above  procedure  to 
TeO?.  Since  the  tellurium  atoms  are  mu:h  heavier  than  the  oxygen  atoms, 
assuming  that  the  tellurium  atoms  form  a  "rigid"  lattice  is  probably  a 
reasonable  approximation.  Using  the  details  from  the  crystal  strucfural 
determination  consider  an  oxygen  with  an  equilibrium  vector  position  (in 

coordinate  units  ot  A)  of  r^  =  (0.63,  1.29,  1.39).  Its  "first"  nearest  neighbor 
tellurium  atoms  are  located  at  r^  =  (3.10,  0.10,  0.00)  and  =  (2.30,  2.49, 
1.91),  which  correspond  to  d^,-,  =  1.92  1  and  d2f|  =  2.09  respectively.  Tne 
bond  angle  between  the  oxygen  and  these  two  tellurium  atoms  is  140. 3°,  which 
differs  considerably  with  the  169°  suggested  in  Reference  22. 

Coupling  the  oxygen  atom  to  only  these  two  tellurium  atoms  is  not 
sufficient  to  constrain  the  oxygen  to  the  specified  equilibrium  position.  At 
least  one  additional  coupling  must  be  included  for  the  equilibrium  position  to 
be  uniquely  determined.  The  next  nearest  neighbors  to  oxygen  atom  #0  are  two 
pairs  of  oxygen  atoms.  One  pair  has  an  interatomic  distance  of  2.65  A  and  the 
other  pair  is  2.75  A  away.  However,  utilizing  any  of  these  four  oxygen-oxygen 
couplings  would  complicate  Eq.  (94);  that  is  Eq.  (94)  would  become  a  coupled 
differential  equation  since  the  motion  of  one  oxygen  atom  would  now  depend  on 
the  variable  positions  of  other  oxygen  atoms. 

In  order  to  maintain  the  uncoupled  form  of  Eq.  (94),  oxygen  atom  #0 
must  be  viewed  as  interacting  only  with  tellurium  atoms.  Therefore,  the  third 
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coupling  is  made  to  the  next  nearest  tellurium  atom,  located  at  r^  =  (-0.10, 
-0.10,  3.31)  with  djQ  =  2.90  A.  The  spring  between  this  tellurium  atom  and 
oxygen  atom  #0  might  be  considered  as  a  "collective"  spring,  representing  the 
coupling  between  oxygen  atom  # 0  and  all  other  neighbors  except  tellurium  atom, 

*1  and  #2. 

The  infrared  frequencies  of  vibration  in  Te02  have  been  determined  from 
reflectivity  and  Raman  spectra. ^  The  two  highest  energy  modes  should  corre¬ 
spond  to  the  longitudinal  and  transverse  motions  of  oxygen  #0  with  respect  to 
telluriums  si  and  *2.  The  experiment^  has  a  pair  of  high  energy  modes  for  both 
the  longitudinal  and  transverse  cases.  The  longitudinal  frequencies  are  812  and 
720  cm”*  whereas  the  transverse  frequencies  are  769  and  643  cm'*.  The  simpli¬ 
fied  model  assumed  above  would  have  these  pairs  of  frequencies  appear  degenerate. 
Therefore,  for  the  purpose  of  determining  the  spring  constants  in  Te0?,  the 
geometric  mean  of  these  pairs  of  frequencies  is  used,  giving  longitudinal  and 
transverse  modes  which  correspond  to  764.6  cm-1  =  13. Ip  and  703.2  cm'1  =  14. 2y, 
respectively.  Note  that  the  geometric  mean  is  employed  instead  of  the  arith¬ 
metic  average  since  the  former  can  be  done  in  frequency  or  wavelength  and  give 
consistent  results,  whereas  the  latter  can  not. 

In  order  to  simplify  the  fitting  procedure  of  the  spring  constants,  the 
two  springs  between  oxygen  atom  #0  and  tellurium  atoms  *1  and  #2  are  assumed  to 
be  characterized  by  the  same  spring  constant,  even  though  the  interatomic  dis¬ 
tances  are  slightly  different.  Then,  by  obtaining  the  eigenvalues  corresponding 
to  the  matrix  specified  by  Eq.  (101)  for  any  given  set  of  spring  constants,  a 
reasonable  fit  to  the  experimental  absorption  frequencies  was  accomplished.  The 
set  of  spring  constants,  Ki  =  K2  =  197  9/sec2  and  K3  =  270  g/sec^,  yield 
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absorptions  at  761.6  cm"1  =  13. 1 3p ,  699.3  cm-1  =  14.30p  and  169.9  cm"*  =  58.86,,. 
The  highest  energy  vibration  ( 13 . 13u)  is  identified  as  a  longitudinal  mode  since 
the  corresponding  eigenvector  is  only  6.6°  away  from  being  parallel  to  the 
interatomic  distance  vector  between  tellurium  atoms  #1  and  #2.  The  other  two 
modes  ( 14 . 30,i  and  58.86^)  are  transverse  because  their  corresponding  eigen¬ 
vectors  are  nearly  perpendicular  to  this  tellurium  interatomic  distance  vector, 
deviating  by  6.6°  and  0.6°  respectively.  The  lowest  energy  absorption  (53.86.) 
can  also  be  compared  with  the  experiment.--*  There  does  exist  a  transverse  mode 
at  174  cm"*  =  57.5.,  which  agrees  rather  well  with  the  calculated  value. 

The  fina1  step  of  determining  the  displacement  of  oxygen  atom  #0  in 
response  to  a  "dc"  electric  field  is  complicated  by  the  fact  that  the  transverse 
effective  charge  in  Eq.  (102)  is  dependent  upon  the  direction  of  displacement  in 
addition  to  the  crystal  structure.  This  dependence  has  been  evaluated  for  a  few 
specific  crystal  structures ,  ■*  and  a  general  solution  will  be  obtained  in  the 
follow-on  to  this  contract.  Until  the  transverse  effective  charge  is  deter¬ 
mined,  the  ionic  displacive  contribution  to  the  second-order  susceptibility  can¬ 
not  be  calculated  in  Te0?.  Furthermore,  the  sensitivity  of  the  displacement 
vector  with  respect  to  the  "collective"  spring  constant  needs  investigation. 

But  this  too  requires  an  evaluation  of  the  transverse  effective  charge  and, 
hence,  will  be  postponed  to  a  later  date. 

One  last  comment  can  be  made  regarding  the  applicability  of  this  lat¬ 
tice  dynamics  model.  The  relative  motion  in  zincblende  compounds  suggested  by 
Fig.  1  should  be  adequately  described  by  Eqs.  (93)  -  (102).  Consequently,  an 
alternate  expression  to  Eq.  (91)  may  be  obtained  which  depends  directly  on  the 
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infrared  absorption  frequencies  rather  than  the  force  constants,  C0  and  Cj. 
Presumably,  these  experimental  absorption  frequencies  are  more  readily  available 
than  the  latter  force  constants  for  the  zincblende  compounds.  Additionally,  the 
similar  ionic  displacive  motion  that  occurs  in  the  chalcopyri tes  should  also  be 
accommodated  by  Eqs.  (93)  -  (192)  by  permitting  distinct  springs  between  a-c  and 
b-c  atom  pairs. 

2 . 1 . 4 . 5  Symmetry 

The  symmetries  of  the  calculated  first-  and  second-order  susceptibility 
tensors,  determined  by  the  summations  over  bonds  in  Eqs.  (62),  (66)  and  (70)  of 
various  combinations  of  interatomic  distance  vectors,  must  reflect  the  symmetry  of 
the  lattice.  This  point  was  illustrated  for  the  zincblende  structure  in  Eqs. 

(75),  (76),  (86),  (89)  and  (90).  The  symmetry  requirements  on  the  second-order 
susceptibility  are  specified  by  the  point  group  symmetry.7  And  as  mentioned  pre¬ 
viously,  thermodynamics  requires  the  first-order  susceptibility  to  be  symmetric. 
The  formulae  associated  with  the  BOM  must  fulfill  these  symmetry  requirements. 

Examining  the  first-order  susceptibility  expression,  Eq.  (62),  it  is 
obvious  that  the  crystal  geometric  factor  for  each  bond,  d,dj,  is  symmetric  in  i 
and  j.  Therefore,  the  summation  over  bonds  is  also  symmetric  in  i  and  j.  The 
electronic  contribution  to  the  second-order  susceptibil ity,  Eq.  (66),  is  almost 
as  simple.  The  crystal  geometric  factor  for  each  bond  can  be  written  as  djdjdk, 
which  is  symmetric  in  i,  j  and  k,  as  is  the  summation  over  bonds.  However, 
Kleinman's  symmetry5*  implies  that  the  second-order  susceptibility  coefficients, 
arising  from  electronic  processes  where  there  is  no  absorption  or  dispersion, 
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are  symmetric  in  all  three  indices.  Consequently,  the  BOM  electronic 
contribution  to  the  second-order  susceptibility  is  seen  to  be  consistent  wito 
Kleinman's  symmetry  condition.  Finally,  no  general  statement  can  be  made,  at 
this  time,  concerning  the  symmetry  of  the  crystal  geometric  factors  in  Eq.  (70), 
which  describe  the  lattice  contribution  to  the  second-order  susceptibility. 
Presumably,  as  indicated  in  Eq.  (28) ,  these  factors  will  turn  out  to  be  sym¬ 
metric  in  i  and  k  only. 

2.2  Experimental 

The  experimental  efforts  in  this  program  have  focussed  on  two  materi¬ 
als:  Te02  (tellurium  dioxide  or  paratel  1  uri te)  and  T^AsSej  (thallium  arsenic 
selenide  or  IAS).  In  much  the  sane  manner  as  was  shown  in  the  theory,  we  have 
chosen  materials  in  which  w.>  hope  to  be  able  to  experimentally  separate  the 
electronic  and  ionic  contributions  to  the  E-0  effect.  TAS  is  a  ternary 
chalcogenide  semiconductor  with  a  moderate  bandgap  of  -  leV,  Although  there  is 
some  evidence  of  soft  mode  behavior2^  in  this  system,  it  is  thought  that  the 
dominant  mechanism  in  TAS  is  the  electronic  response.  Te02  possesses  character¬ 
istics  that  suggest  comparable  contributions  from  both  effects  to  the  E-0  re¬ 
sponse.  Being  a  well  known  acousto-optic  material ,26  Te02  exhibits  a  soft  mode 
phase  transition26  under  moderate  pressure  (-  9  K  bars).  Additionally  Te02  has 
a  reasonably  high  c(*>)  of  approximately  4  and  an  e(o)  of  about  25.  Our  experi¬ 
mental  efforts  in  TljAsSe^  to  date  have  been  of  a  rather  preliminary  nature  and 
will  be  completed  during  the  next  few  months.  The  Te02  work  has  progressed  much 
more  quickly  and  is  essentially  complete  (see  Appendix  A)  at  present. 
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2.2.1  TeO? 

The  point  group  symmetry 22  of  TeC^  is  tetragonal:  422.  Tnis  implies 
that  the  material  is  optically  uniaxial  and  has  only  two  non-zero  electro-optic 
coefficients  (r^j  =  -r^).^  TeO^  is  also  weakly  piezoelectric  having  one  non¬ 
zero  piezoelectric  coefficient^  d^-  =  8.13  x  10-1^  C/N,  indicating  there  is 
considerable  response  of  the  lattice  to  the  presence  of  a  dc  electric  field. 
Although  Kleinnan's  symmetry^  suggest  that  materials  having  the  symmetry  of  422 
should  have  non-zero  second-order  susceptibil ity  coeff icients ,  the  underlying 
assumption  of  Kleinnan's  formalism  make  it  inappropn ate  when  one  of  the  elec¬ 
tric  fields  is  dc  (j  =  0),  as  is  pointed  out  above  (see  Theory).  In  any  case, 
experiments14, 29  yiave  s^own  that  there  is  an  apparent  "violation"  of 
Kleinnan’s  symmetry  in  Te02  even  at  optical  frequencies  (both  electric  fields 
have  u  *  0).  We  shall,  therefore,  examine  the  details  of  the  dc  E-0  effect  in 
this  material . 

If  light  is  propagated  along  the  <01 1 >  direction,  the  two  normal  modes 
occur  with  the  polarization  along  <100>  (ordinary)  and  along  <0ll>  (extraordi¬ 
nary).  The  impermeability  tensor  becomes: 


/  1 


B  = 


o 

0 


-  r,,E  \ 

+>  y 


+  r E 
41  x 


\  -r4i‘ 


*  r41Ex 


VI 


If  we  write  the  polarization  vector: 


(103) 


P  = 


(104) 
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where  =  the  polarization  direction  cosines,  we  can  calculate  the  dielectric 
response  by: 


-  R  =  P(B)P  = 
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Now  in  particuid'  for  E  =  E<5  and  P  =  (100),  which  is  the  ordinary 
pol arization: 
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where  ng(eff)  is  the  effective  extraordi nary  index  for  a  propagation  direction 
of  <0 1 1 > .  Thus  it  is  only  the  extraordinary  index  (or  mode)  that  is  affected 
for  our  particular  choice  of  the  crystal  cut  which  is  shown  in  Fig.  3. 

An  experimental  sample  which  utilized  the  interaction  geometry  shown  in 
Fig.  3  was  constructed  from  Tef^  single  crystal  material  that  was  obtained 
commercially  from  Crystal  Technology  in  Mountain  View,  California. 

The  actual  dimensions  were  1.9  x  0.57  x  0.72  cm.  Crystallographic 
alignment  was  accomplished  using  the  X-ray  back  reflection  Laue  technique.  A 
typical  alignment  shot  taken  normal  to  the  <100>  face  is  shown  in  Fig.  4.  The 
sample  was  wire  sawed  out  of  a  larger  boule  and  then  ground  and  polished.  After 
final  polishing  of  the  end  optical  faces,  the  sample  was  viewed  using  monochro¬ 
matic  (632BA)  light,  as  shown  in  Fig.  5.  When  monochromatic  polarized  light  was 
propagated  to  excite  both  optical  modes,  and  then  examined  at  the  output  with  a 
crossed  polarizer,  an  interference  pattern  (see  Fig.  6)  resulted  which  indicated 
that  the  sample  was  slightly  wedged.  By  simultaneously  reflecting  a  laser  beam 
off  tne  front  and  back  surfaces,  the  wedge  angle  was  found  to  be  approximately 
0.3  degrees.  This  agrees  well  with  the  number  that  is  derived  by  measuring  the 
distance  between  the  fringes,  and  substituting  into  the  following  relations.  The 
phase  difference,  r,  between  the  ordinary  and  extraordinary  mode  is  given  by: 

r  = -~  (ne(eff)  -  nQ)  (108) 

where  X  =  optical  wavelength,  ne(eff)  and  nQ  are  the  refractive  indices,  and  a  = 
the  interaction  length.  The  phase  difference  between  adjacent  fringes  is  given 
by: 
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Fig.  3  TeOo  electro-optic  sample  configuration 
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Ar  =  (ne(eff)  -  nQ) 


(109) 


where  and  z ^  are  the  respective  interaction  lengths  for  the  two  fringes.  Tne 
absolute  phase  difference,  between  adjacent  light  (or  dark)  fringes  a  distance 
dg  apart,  is  rt.  Further,  one  can  relate  the  distance  between  fringes,  dg,  to 
the  wedge  angle,  a,  by: 


do=  (un) 

using  the  smal 1  angle  approximation.  If  the  optical  path  length  changes,  the 
fringe  pattern  will  move  up  or  down  depending  on  the  sign  of  the  effect.  Of 
course,  this  change  in  optical  path  length  may  result  from  either  of  two  different 
effects:  1)  a  change  in  the  physical  length  of  the  crystal  due  to  mechanical 
strains  induced  by  the  inverse  piezoelectric  effect  and  2)  a  change  in  the  (ex¬ 
traordinary)  index  of  refraction.  We  shall  show  below  that  the  first  effect  is 
negligible  in  Tet^.  Therefore,  assuming  that  the  fringes  move  only  as  a  result 
of  a  change  in  the  index,  the  distance  the  fringe  pattern  moves,  dj,  is  given 
by: 


„  SL  (eff)  -  ne(eff)) 

V  5  ~(v7erfT~%T—  011 

in  the  limit  of  small  a.  Here  n  '(eff)  =  n  (eff)  +  An  (eff)  with  An  (eff)  = 
-2ngJ(eff)  r^Ex  where  we  have  assumed  that  the  index  has  been  modified  by  a  dc 
electric  field,  Ex.  Taking  the  ratio  of  d^/dg  and  solving  for  r^,  we  obtain: 
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rJL  =  A(d1/d0)/(4Exne3(eff) 


(112; 


All  that  is  required  then  is  to  measure  the  fractional  shift  in  the  f-inge 
pattern  for  a  given  applied  electric  field.  In  order  to  quantitatively  deter¬ 
mine  this  ratio,  the  sample  was  placed  in  the  optical  setup  shown  in  Fig.  5. 

Photographs  of  the  fringe  patterns  and  the  fiducial  crosshairs  for 
various  applied  fields  are  shown  in  Figs.  6a  -  6c.  The  applied  voltages  were 
-4.9,  0  and  *  4.9  kV.  The  result  was  d^/dQ  »  1/2  for  9.3  k V.  This  corresponds 
to  an  index  change  of  pne(eff)  =  8.3  x  10"^  and  an  electro-optic  coefficient  of 
lr4l!  =  0.76  x  10'12  m/V.  The  sign  of  this  electro-optic  coefficient  was  ob¬ 
tained  by  correlating  the  direction  of  the  fringe  shift,  the  wedge  angle  and  the 
sign  of  the  applied  voltage.  The  result  is  that  r4l  =  -  0.76  x  lO-1^  m/V.  The 
closely  spaced  fringes  in  the  lower  right  corners  of  Figs.  6a  -  6c  are  probably 
due  to  a  slight  rounding  of  the  presumed  flat  optical  face  near  the  edge  of  the 
sample. 

An  estimate  of  the  physical  distortion  to  the  sample  due  to  the  strain 
resulting  from  the  inverse  piezoelectric  effect,  is  given  by  the  following: 


e .  =  \  d  . .  E  . 
1  j  J1  J 


(113) 


where  is  the  strain,  is  the  piezoelectric  coefficient  and  Ej  is  the  ap¬ 
plied  electric  field.  The  only  non-zero  piezoelectric  coefficient00  for  the 
point  group  422  is  d14.  The  measured  value27  for  Te02  is  d14  =  8.13  x  10*12  C/N. 
The  change  in  wedge  angle  resulting  from  the  piezoelectric  strain  induced  by  an 
electric  field  of  9.8  kV  is  0.009°  or  0.3%  and  will,  therefore,  be  neglected. 
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The  strain  induced  by  the  inverse  piezoelectric  effect  does,  however, 
affect  the  index  via  the  strain-optic  effect  giving  rise  to  a  change  in 
impermeabi  1  ity  given  by:  31 


*Bi  =  l  PijeJ  =  AVkA 

J  J  ♦  * 


(114) 


where  P^j  is  the  relevant  component  of  the  photoelastic  tensor  for  Tef^.  Indeed 
we  may  define  a  secondary  electro-optic  coefficient  by: 


P,  .d. 
ij 


(115) 


Since  in  our  experiment  EK  =  Ej  and  d^j  =  d^,  P44  is  the  only  non-zero 
component  which  contributes  to  the  secondary  E-0  effect.  This  component  has 
been  measured^  and  has  the  value  P44  =  -0.17.  The  secondary  E-0  coefficient 
therefore  has  the  value: 


r4f  =  p44di4  =  “1-38  *  10'12  m/V  (116) 

Using  the  applied  dc  voltage  of  9.8  kV,  we  deduce  a  change  in  effective  index  of 
An  (eff)  =  +1,50  x  10'5  .  (117) 


If  we  now  subtract  this  result  from  the  number  which  was  experimentally  derived 
(at  constant  stress),  we  obtain 
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(An  (eff)) 


exp 


[Ane(eff))s!r.0  =  8.3  x 


10'°  -  1.5  x  10 


-5 


=  -0.67  x  10"5 


(118) 


This  in  turn  yields  a  value  for  the  primary  E-Q  coefficient  at  constant  strain 
of : 


r 


S 

14 


+0.62  x 


m  /  V 


(118) 


where  the  superscript  S  indicates  constant  strain  ("clamped").  This  value  can 
be  compared  with  other  nonlinear  optical  measurements  in  Te02;  for  example, 


dGJGj  =  1.45  x  10"9  esu  =  iox  0.61  x  10‘12  m/V 


(120) 


derived  from  second  harmonic  generation  measurements.  14,15,28,29  (i\|0te  that  the 

d^  here  is  distinct  from  the  piezoelectric  coefficient  referred  to  above.) 

These  comparisons  should  not  be  made  directly  between  r^  and  d^J.  Instead  they 
should  be  made  between  the  electro-optic  susceptibility: 


eo 

*41 


n2n2r  $ 
o  e  41  _  .  . 

TtT - 1,5 


10"12  m/V 


(121) 


and  the  second  harmonic  generation  susceptibility  d 


SHG 
14  • 
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2.2.2  Tl^AsSe, 

Thallium  arsenic  selenide  or  TAS  is  a  ternary  semiconductor  which 
crystallizes  with  the  trigonal  point  group  symmetry  3m.  This  implies  that  TAS 
is  optically  uniaxial  and  should  have  4  non-zero  electro-optic  components:7  r,^, 
r22’  r33 »  arK*  r5i*  This  compound  has  been  the  subject  of  extensive  investiga¬ 
tion  for  its  optic, ^  acousto-optic,^"^  and  non-linear^  (second  harmonic 
generation)  properties.  It  is  transparent  from  the  near  (0.9  ym)  to  mid  (14  um) 
IR.  By  using  some  of  the  formalism  associated  with  the  dielectric  impermeabil¬ 
ity  tensor  derived  above  (Section  2.2.1),  we  now  examine  the  electro-optic 
effect  in  TAS. 


As  a  preliminary  example,  let  us  consider  the  case  of  an  x-y-z 
rectangular  prism  cut  so  that  the  long  dimension  is  along  x.  The  dc  electric 
field  is  applied  normal  to  the  long  direction  and  along  y.  The  third  direction 
or  z,  is  parallel  to  the  optic  axis,  light  is  then  propagated  along  x  and 
polarized  along  y  (ordinary)  or  along  z  (extraordinary).  The  dielectric 
impermeability  tensor  is  then: 


B  = 


1 

o 

0 


+  r22Ey 
r51Ey 


(122) 


The  off-diagonal  components,  r^Ey,  produce  a  small  rotation  of  the 
indicatrix  about  the  x  axis  by  an  amount: 


tan  2o 


*51^ 


1  K  ♦  r 

"o  "e 


22  Ey 


(123) 
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An  upper  estimate  of  this  rotation  (assuming  r^j  =  100  x  10"^  m/V)  for  voltages 
of  5KV  is  9  <  1°,  and,  therefore,  this  effect  will  be  ignored.  Light  polarized 
along  y  (ordinary)  has  the  effective  index: 


n 


nor22Ey 

22 


(124) 


while  the  other  polarization  (along  z  or  extraordinary)  remains 


n  =  ne  .  (125) 

Therefore,  it  is  only  the  ordinary  component  which  is  affected  by  the  dc 
field.  If  we  then  arrange  to  have  the  incident  light  linearly  polarized  45° 
between  y  and  z,  or  along  <01 1>,  both  modes  will  be  excited.  Due  to  the 
difference  in  propagation  velocities,  the  two  modes  will  emerge  out  of  phase 
with  one  another  by  an  amount,  r,  given  by: 

r  -  M  (An)  (126) 

where  as  before  %  is  the  interaction  length,  x  is  the  wavelength  and  in  is  the 
birefringence  (ne-nQ).  The  polarization  state  at  the  exit  face  will  be  in 
general  ellipitcal,  with  circular  and  linear  polarizations  as  specific  cases 
(for  90°  and  0°  or  180°  phase  differences,  respectively).  As  we  have  shown,  if 
an  electric  field  is  applied  along  y,  then  the  birefringence,  and  hence,  r, 
changes  by: 
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(127) 


The  total  phase  retardation  is  then 


+  AT 


-If  ((An; 


-  n^E  r  ,] 
o  y  22' 


(12Hi 


To  experimental ly  observe  this  effect,  we  set  up  the  sample  and  optics 
as  shown  in  Fig.  7.  The  laser  used  was  a  He-Ne  3.39  yn  laser.  The  compensator 
was  adjusted  to  give  linearly  polarized  light  with  no  field  applied.  The  pola¬ 
rizer  was  then  "crossed"  so  that  essentially  na  light  was  observed  at  the  detec¬ 
tor.  As  the  electric  field  was  increased,  a  series  of  oscillations  in  the  sig¬ 
nal  from  the  detector  was  observed  as  the  state  of  the  light  exiting  the  crystal 
changed  ellipticity.  The  voltage  required  to  change  the  signal  by  one  oscilla¬ 
tion  (or  the  phase  by  k)  is  called  the  half-wave  voltage  Vn  and  is  given  by: 


AF 


n’Er,, 
o  22 


(129) 


TT 


V  r 

TT  CL 


— -  -Ty— 


(130) 


where  y  is  the  distance  between  the  electrodes.  Letting  y  and  i  be  of  unit 
dimension : 


V 

u 


\ 

n7r 


o'  22 


(131) 
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The  sample  used  in  this  preliminary  experiment  was  fashioned  from  an 
in-house  grown  » 37  boule  of  Tl3AsSe3«  The  sample  was  first  wire  sawed  from  the 
boule  and  then  mechanically  ground  and  polished  (on  the  optical  faces)  to  fin¬ 
ished  dimensions  of  0.36  x  0.32  x  1.27  cm  for  y,  z  and  x,  respectively.  Silver 
paint  was  used  for  electrodes.  The  sample  resistance  with  no  applied  field  was 
greater  than  20  Mohm  which  meant  the  resistivity  was  more  than  2  x  107  ohm-cm. 
This  is  consistent  with  the  room  temperature,  and  an  approximately  1  eV  band- 
gap.36  Figure  8  shows  the  data  for  this  preliminary  experiment. 

Using  the  voltage  between  peaks  of  2.2  kV  and  the  sample  dimensions  of 
1.27  cm  for  z  and  0.36  cm  for  y  and  the  room  temperature  index  at  3.39  ;)n  of 
3.37,  we  find  r^  =  11  x  10"^  m/V,  and  the  half  wave  voltage  V ^  =  8  k V.  This 
result  includes  the  piezoelectric-induced  effects. 

Our  next  experiment,  to  be  carried  out  under  the  follow-on  to  this 
contract,  will  be  to  place  the  electric  field  along  z,  with  the  same  sample 
geometry  as  before.  The  dielectric  impermeability  tensor  then  becomes 


B  = 


n?  +  r13E3 
o 

0 

0 


h  +  r  13  E  3 


W  +  r33E3  / 


(132) 


With  light  propagating  along  x  and  polarized  along  y,  (the  ordinary  mode),  the 
index  becomes: 


n'  =  n 
o  o 


0  13  3 


(133) 
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The  orthogonal  polarization  along  z  (the  extraordinary  mode)  has  the  followin' 
i ndex  : 


O' 

e 


n^r^E., 
e  33  3 


(13i; 


It  is  common  practice  to  define  the  effective  birefringence  in  terms  of  a 
combined  electro-optic  coefficient  rc: 


with 


(n  -  n  )  - 
'  e  o' 


"7r —  n 3r -  n - r , . 
2  K  e  33  o  13 


E3rcn| 

(n  -  n  )  - x — 

'  e  o'  2 


r  1  >  "o'*13 

c  '  2  33  '  ~"nj 

e 


The  field  dependent  phase  retardation  is  then 


X  1  \  ner33 


nor13)E: 


~  i  r  E,n3 
X  c  3  e 


(135) 


(136) 


(137) 


The  half-wave  voltage  in  this  case  is: 

2"  ■  T  1  rcVl/» 


(138) 
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V 

IT 


n'r 
e  c 


(139) 


for  unit  sample  dimensions. 

It  will,  therefore,  be  possible  to  measure  rc  using  the  same  experi¬ 
mental  arrangement  as  before,  by  merely  removing  the  silver  paint  electrodes 
from  the  y-faces  and  placing  them  on  the  z-faces. 

The  coefficient  r^  can  be  independently  obtained  by  a  longitudinal 
measurement  on  a  z-cut  plate.  The  index  change  seen  by  light  propagating  along 
the  z  direction  polarized  in  the  x-y  plane  (pure  ordinary)  is  simply: 


n' 


n  1 

no  ■  -r  r13E3 


(140) 


which  can  be  measured  in  a  Fabry-Perot  configuration.  Then  by  knowing  r13  and 
rc,  we  can  calculate  r33.  The  coefficient  r^  can  be  measured  as  a  rotation  of 
the  indicatrix  using  a  longitudinal  geometry  with  propagation  along  the  y  axis. 

Finally,  in  order  to  obtain  the  "clamped"  electro-optic  coefficients, 

C 

r?  ,  it  will  be  necessary  to  experimentally  determine  the  piezoelectric 
coefficients  and  photoelastic  constants  of  Tl3AsSe3. 
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3.0  SUMMARY  AND  RECOMMENDATIONS 

In  summary,  this  program  has  advanced  well,  with  significant  progress 
being  made  on  both  the  theoretical  and  experimental  portions  of  this  program. 

The  groundwork  has  been  established  for  the  theoretical  approach  to  the  electro¬ 
optic  effect  in  solids.  The  underlying  philosophy  to  our  method  is  to  use  a 
generalized  "chemical  bond"  approach  to  the  first-  and  second-order  electric 
susceptibility.  This  methodology  utilizes  nearest-nei ghbor  interactions,  an! 
parametrizes  the  relevant  matrix  elements  in  terms  of  universal  atomic  param¬ 
eters.  Calculations,  therefore,  only  require  a  knowledge  of  elemental 
composition  and  the  atomic  bond  configurations.  This  approach  lends  itself  to 
predictive  comparisons  of  electro-optic  response  not  only  within  a  given  crystal 
structure  (making  elemental  substitutions)  but  between  different  structures  as 
well. 

Details  of  the  methodology  associated  with  the  Rond  Orbital  Model  have 
been  given  along  with  applications  to  the  specific  cases  of  the  chalcopyrite  and 
zincblende  crystal  structures.  Initially,  however,  the  Bond  Orbital  Model  was 
only  applicable  to  calculations  of  the  response  of  the  electronic  system  to 
probing  fields  of  frequencies  above  the  Restrahl  frequency  and  below  the  first 
inter-band  or  interactomic  transition.  Therefore,  the  next  level  of  sophis¬ 
tication  introduced  into  the  problem  was  an  attempt  to  include  the  lattice 
dynamics. 

The  approach  here  (still  being  developed)  was  formulated  in  terms  of 
the  microscopic  spring  constants.  Knowing  these  spring  constants,  one  could,  in 
principle,  calculate  the  relative  displacement  and  resulting  lattice  distortion 
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of  a  crystal  in  the  presence  of  a  dc  electric  field.  Then,  by  knowing  the  new 
equilibrium  ionic  positions,  one  can  recalculate  the  susceptibility  and  hence 
the  electro-optic  effect.  The  process  is  evolving  along  lines  oriented  towards 
obtaining  these  microscopic  spring  constants  from  macroscopic  measurable  proper¬ 
ties  of  the  materials  such  as  elastic  constants  or  lattice  mode  frequencies. 

This  approach  tends  to  somewhat  diminish  the  overall  generality  and,  therefore, 
predictive  properties  of  the  method  in  that  more  measurements  are  required  be¬ 
fore  the  calculations  may  be  undertaken.  This,  however,  should  not  be  taken  as 
an  end  result,  but  rather  as  an  intermediate  step  in  an  evolving  process.  At 
some  future  date,  it  is  hoped  al so  to  be  able  to  generalize  the  lattice  response 
to  structure  and  elemental  composition,  thereby  recovering  the  predictive 
capabilities. 

The  experimental  portion  of  this  program  has  involved  measurements  of 

the  electro-optic  effect  in  TeO^  and  Tl^AsSej.  The  TeO^  data  is  complete  and 

S  - 12 

has  yielded  thr>  val  je  r^  =  -  9.61  x  10  m/V.  Experiments  on  TljAsSej  at 

T 

3.39  um  ha/e  just  begun  but  a  preliminary  result  indicates  jr^l  =  11  x  10  m/V. 
Further  measurements  will  be  carried  out  on  Tl^AsSe-j  during  the  next  phase  of 
the  program  to  obtain  and  r^.  Additionally,  other  materials  such  as  SbSl 
and  Barium  Titanium  Niobate  (BTN)  will  be  approached  experimental ly  during  the 
next  phase  of  this  program.  SbSI  has  one  of  the  highest  reported  piezoelectric 
coupling  coefficients  of  any  solid,*®  and  should  have  a  correspondingly  large 
electro-optic  effect.  BTN  has  been  shown  to  have  a  significant  electro-optic 
effect.-*9  Structural  considerations  have  shown  that  it  may  be  possible  to 
substitute  lanthanum  or  potassium  on  some  of  the  sites  in  the  unit  cell. 


66 

C3485A/jbs 


Rockwell  International 

Science  Center 


SC5266.1FK 

Our  recommendation,  therefore,  is  to  proceed  with  the  theoretical  w on¬ 
to  approach  a  completely  generalized  formalism  which  has  predictive  capabili¬ 
ties.  In  addition,  the  experimental  work  should  proceed  to  evaluate  materials 
which  provide  insight  and  feedback  to  the  theory,  as  well  as  provide  new  high 
performance  materials  for  further  applied  device  efforts. 
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APPENDIX  I 


LINEAR  ELECTRO-OPTIC  EFFECT  IN  TELLURIUM  DIOXIDE* 

M.D.  Ewbank  and  P.R.  Newman 

Rockwell  International  Science  Center 
1049  Camino  dos  Rios 
Thousand  Oaks,  CA  91360 


ABSTRACT 


The  only  non-zero  dc  electro-optic  (or  Pockel ’ s)  coefficient  in 
tellurium  dioxide  was  measured  at  constant  stress.  The  result  was  r^  =  -0.76 
x  10'12  m/V.  The  secondary  electro-optic  effect,  due  to  the  converse 
piezoelectric  and  photoelastic  effects,  is  approximately  twice  the  observed 
electro-optic  response. 


♦This  work  was  supported  by  ONR  contract  number  N00014-80-C-0498. 


1 

J3376A/bw 


I.  INTRODUCTION 


Tellurium  dioxide  (paratel lurite  or  Te02)  is  a  well-known  acousto- 
optic  material,  utilized  in  devices  such  as  acousto-optic  tunable  fil¬ 
ters,^’4  beam  deflectors  and  modulators.  The  anomalously  slow  acoustic  shear 
velocity, *  in  the  <110>  direction  which  can  be  associated  with  the  onset  of  a 
soft  lattice  mode  transition,^  leads  to  a  favorable  acousto-optic  figure-of- 
merit.  The  presence  of  nearby  phase  transitions  in  ferroelectric  materials 
which  exhibit  superior  electro-optic  behavior^  suggests  that  the  electro-optic 
effect  in  TeC^  may  also  be  enhanced  by  this  softening  lattice  mode.  Conse¬ 
quently,  the  room  temperature  value  of  the  single  non-zero  dc  electro-optic 
coefficient  (r^)  for  TeC^  was  determined  experimentally  at  a  wavelength  of 
6323  A. 


II.  EXPERIMENT 

The  point  group  symmetry7  of  tellurium  dioxide  is  422,  which  implies 
that  Te02  is  optically  uniaxial  and  has  only  two  non-zero  electro-optic 

Q 

coefficients  (r^  =  -rg^).  An  optical  quality  sample  of  TeC^  was 
fabricated,  using  material  obtained  commercially  from  Crystal  T~-  olo  ,  Inc-,, 
with  the  orientation  for  the  measurement  of  r^  indicated  in  Fig  1 .  -> 

optical  faces,  1.9  cm  apart,  were  perpendicular  to  the  <01 1 >  direction,  and 
the  silver-painted  electrodes,  separated  by  0.57  cm,  were  normal  to  <100>.  A 
schematic  of  the  optical  setup  is  illustrated  in  Fig.  2.  At  a  wavelength  of 
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6328  A,  Te02  is  optically  transparent  and  its  refractive  indices  at  room 
temperature^  are  n0  =  2.2585  and  ng  =  2.4112. 

For  light  propagating  in  the  <01 1 >  direction  and  polarized  in  the 
<2ll>  direction  (i.e.,  at  45  degrees  fron  <100>  and  <0 I 1 > )  both  the  ordinary 
and  extraordinary  modes  of  propagation  are  equally  excited.  When  the  ordinary 
and  extraordinary  beams  exit  the  Te02  sample,  the  relative  phase  shift,  r, 
between  ordinary  and  extraordinary  beams,  assuming  that  the  dc  electric  field 
is  zero  and  that  the  front  and  back  surfaces  are  exactly  flat  and  parallel, 
will  be 


r  =  (2iU/A)(ne(eff)  -  nQ) 


(1) 


where  z  is  the  sample  interaction  length,  A  is  the  wavelength,  and  n0  is  the 
ordinary  index.  The  effective  index  of  the  extraordi nary  wave,  ng(eff),  for 


this  propagation  direction  is  given  by  the  relation 


10 


ne(eff)  =  n  n0ne//nj  +  V 


(2) 


with  ng  being  the  extraordi nary  index. 

Since  the  ordinary  and  extraordinary  waves  have  equal  amplitudes,  the 
intensity  of  the  analyzed  beam  will  vary  with  the  orientation  of  the  analyzing 
polarizer,  because  the  polarization  state  of  the  unanalyzed  light  is  either 
linear,  elliptic  or  circular  depending  on  the  relative  phase  retardation  of 
the  two  modes.  If  the  front  and  rear  surfaces  of  the  sample  are  wedged 
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slightly,  an  interference  fringe  pattern  will  result.  The  spacing,  dg, 
between  adjacent  fringes  which  have  a  relative  phase  shift  of  can  be 
written: 


do  5  Mne(eff)  -  n0J  (3> 

in  the  limit  that  the  wedge  angle,  a,  is  small. 

When  an  external  dc  electric  field,  Ex,  is  applied  in  the  <100> 
direction,  the  extraordinary  wave  experiences  a  change  in  impermeability 
of  r4i ,  which  corresponds  to  an  effective  extraordinary  index  change, 
Ang(eff),  given  by10 

Ane(eff)  s  -rj^n^eff)/^  (4) 

in  the  approximation  that  this  change  in  index  is  much  less  than  the  index 
itself.  The  superscript  "T"  indicates  that  the  measurement  is  performed  under 
constant  stress  ("undamped")®.  (This  change  in  effective  index  is  maximized 
for  the  <0 1 1 >  propagation  direction  only  in  the  limit  of  the  birefringence 
approaching  zero.  For  the  birefringence  of  Te02  at  6328  A,  the  maximum  occurs 
when  propagating  about  3  degrees  away  from  the  <01 1>  direction;  but  the 
difference  in  effective  index  change  for  these  two  propagation  directions  is 
less  than  0.5%  and  therefore  is  negligible.)  Since  the  dc  electric  field  does 
not  affect  the  ordinary  wave,  the  phase  shift  between  ordinary  and 
extraordinary  modes  will  change  with  external  electric  field  by  an  amount  that 
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can  be  determined  by  combining  Eqs.  (1)  and  (4).  More  specifically,  the 
fringe  pattern  will  shift  by  a  distance,  d},  when  the  dc  field  is  applied: 


*l>g(eff)  -  ne(eff}] 
d!  =  “  aLn^(eff)  -  nQJ 


(5) 


where  n^(eff)  is  the  effective  extraordinary  index  for  a  non-zero  dc  field 
(i.e.,  n^eff)  =  ng(eff)  +  Ang(eff)).  The  fringe  shift  with  respect  to  the 
separation  between  fringes  can  then  be  written  as  the  ratio: 

(dx/d0)  s  2*  Ane(eff)/X  (6) 

in  the  approximation  that  the  change  in  effective  extraordinary  index  is  much 
less  than  the  effective  bi refringence,  [ng(eff)  -  nQ].  Finally,  the  electro¬ 
optic  coefficient,  rjp  can  be  expressed  as  a  function  of  this  fractional 
fringe  shift,  the  external  dc  electric  field,  the  wavelength,  the  sample 
interaction  length  and  the  effective  extraordinary  refractive  index  in  the 
following  form: 

t  ~ ^  ( d-j  /  dp, ) 

r41  2  ZExn|(eff5  *  (7) 

Voltages  of  -4.9,  0.0  and  +4.9  kV  were  supplied  to  the  electrodes  on 
the  x-faces  of  the  Te02  sample  and  produced  the  interference  fringe  patterns 
shown  in  Fig.  3.  The  relative  motion  of  the  fringes,  or  fractional  fringe 
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shift,  was  determined,  from  photographic  measurements  referenced  to  a  fixed 
cross-hair,  to  be  approximately  1/2  fringe  (±  20%)  for  an  applied  voltage 
change  of  9.8  kV.  It  should  be  noted  that  the  closely  spaced  fringes  in  the 
lower  right  corners  of  Figs.  3a,  b,  and  c  are  probably  due  to  a  slight  round¬ 
ing  of  the  presumed  flat  optical  face  near  one  edge  of  the  sample.  Addi¬ 
tionally,  the  fine  rectangular  fringe  pattern  is  caused  by  Fresnel  diffraction 
from  the  sample  aperature.  The  wedge  angle  was  found,  by  comparing  simul¬ 
taneous  reflections  from  both  optical  faces,  to  be  approximately  0.3  degrees, 
which  is  roughly  consistent  with  the  number  of  major  fringes  shown  in  Fig.  3. 
Since  the  wedge  angle  is  so  small,  Eq.  (7)  has  been  used  to  calculate  a  room 
temperature  value  for  the  magnitude  of  the  electro-optic  coefficient  of  Te02 
at  6323  A,  which  was  |r^|  =  0.76  x  10'^  m/V.  Equation  (4)  indicates  that 
this  corresponds  to  an  effective  extraordi nary  index  change  of  Ang(eff)  =8.3 

x  10-6. 

The  sign  of  this  electro-optic  coefficient  was  obtained  by  correla¬ 
ting  the  direction  of  fringe  shift  with  respect  to  the  sign  of  applied  voltage 
and  the  orientation  of  the  wedge  angle  with  respect  to  the  optic  axis.  The 
orientation  of  the  optic  axis  was  determined  by  noting  that,  for  this  direc¬ 
tion  of  propagation,  the  extraordinary  beam  experienced  a  "walk-off",  with 
respect  to  the  ordinary  beam,  of  3.7  degrees  in  the  direction  of  the  optic 
axis.  Referring  to  the  fringe  patterns  in  Fig.  3,  the  electrodes  lie  on  the 
faces  perpendicular  to  the  fringes,  with  the  ground  electrode  being  on  the 
upper  right  face.  The  wedging,  as  determined  by  retroreflection,  indicated 
that  the  interaction  length  in  the  portion  of  the  sample,  which  corresponds  to 
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the  region  near  the  upper  left  surface  of  Fig.  3,  is  longer  than  the  interac¬ 
tion  length  near  the  lower  right  surface.  When  a  positive  voltage  was  applied 
to  the  lower  left  face,  the  fringes  shifted  toward  the  lower  right  where  the 
interaction  length  has  decreased.  In  order  to  maintain  a  constant  phase 
difference  for  a  given  fringe,  the  optical  path  length  must  remain  the  same. 
Since  the  physical  path  length  decreased  for  a  positive  voltage,  the  index 
apparently  increased  to  yield  a  constant  optical  path  length.  Then,  by 
utilizing  Eq.  (4),  the  electro-optic  coefficient,  rL,  must  be  negative. 

Finally,  the  results  of  this  electro-optic  measurement  on  TeC^  at 
63^3  A  can  be  summarized  as 

rjj  =  -0.76  x  10'12  n/V  (8) 

for  the  room  temperature  value. 


III.  DISCUSSION 

The  above  experiment  with  TeOj  measures  relative  changes  in  optical 
oath  lengths  between  ordinary  and  extraordinary  waves,  caused  by  the  applica- 
* • on  of  an  external  dc  e  trie  field.  However,  thre"  physical  phenomena  can 
if*  a  change  in  the  optical  path  length:  (1)  the  electro-optic  effect 
„*  ■»  ctly  changes  the  refractive  index,  (2)  the  converse  piezoelectric 
•  .wch  changes  the  physical  dimensions  of  the  crystal  (in  particular, 

.«•  and  (3)  the  secondary  electro-optic  effect  (the  photoelastic 
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or  strain-optic  effect  due  to  the  converse  piezoelectric,  induced  strain) 
which  also  modifies  the  refractive  index. 

The  strain,  e^ ,  which  results  from  the  converse  piezoelectric  effect, 
can  be  calculated  by  using  the  relation:11 


ei 


=  1  d 


ji  J 


(9) 


where  d^  is  the  piezoelectric  coefficient  and  E^  is  the  electric  field.  The 
physical  deformation  can  then  be  estimated  from  the  measured  value  of  the  only 
non-zero  piezoelectric  coefficient1^  (d^  =  8.13  *  10-1^  C/N)  and  the  external 
dc  electric  field.  This  shear  strain  corresponds  to  a  change  in  wedge  angle 
of  approximately  0.3%  for  an  applied  voltage  change  of  9.8  kV  and  a  negligible 
fringe  shift. 


The  secondary  electro-optic  effect  gives  rise  to  a  change  in 

13 

impermeability,  A8^,  given  byA 


A8 


i  '  j  Pij  6j  "  kPi j  dkj 


(10) 


where  p^  is  the  photoelastic  constant,  ej  is  the  piezoelectric  strain,  dkj  is 

the  piezoelectric  coefficient  and  E^  is  the  electric  field.  Then,  the 

(2) 

secondary  electro-optic  coefficient.  r>k',  can  be  defined  by  the  equation 


^  Pij  dkj 


(ID 
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Since  d^4  is  the  only  non-zero  piezoelectric  coefficient,  only  P44  contributes 
to  rjjp  in  Eq.  (11),  and  this  photoelastic  constant  has  been  found 
experimentally1  (p44  =  -0.17).  The  secondary  electro-optic  coefficient  is 
-1.38  x  10"iC  m/V,  which  corresponds  to  a  change  in  effective  extraordinary 
index  of  Ang(eff)  =  +1.50  x  10'^  using  an  equivalent  form  to  Eq.  (4)  with  a 
1.7  x  106  V/n  electric  field.  Subtracting  this  secondary  electro-optic 
coefficient  from  the  measured  value  at  constant  stress  (see  Eq.  (8)),  one 

1  O 

obtains  the  primary1'1  electro-optic  coefficient  for  TeO ?: 

r^  »  +0.62  x  IQ'12  n/V  (12) 

where  the  superscript  "S"  indicates  constant  strain  ("clamped").6 

There  have  been  numerous  measurements  of  the  second  harmonic 
generation  (SHG)  coefficient  in  TeO.?9, 14-16  which  have  shown  a  violation  of 
Kleinman's  symmetry  relation.1^  The  results  of  these  measurements  yield  a 
reasonably  consistent  value  of  approximately 

Id,  J  =  1.45  x  10"9esu  =  e  x  0.61  x  10-12  m/V  (13) 

for  a  fundamental  wavelength  of  1.064  urn,  and  where  eQ  is  the  permittivity  of 
free  space.  Note  that  this  d^  is  distinct  from  the  piezoelectric  coefficient 
in  Eqs.  (9)-(ll).  In  addition,  both  the  symmetric  and  antisymmetric  nonlinear 
optical  susceptibilities  for  sum  and  difference  mixing  have 
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been  measured  in  TeC^,  and  the  symmetric  coefficient  is  in  good  correspondence 
with  the  SHG  value  in  Eq.  (13). ^ 

Previous  attempts  have  correlated  the  second  order  susceptibilities 
of  optical  rectification  and  the  linear  electro-optic  effect.19  Also,  com¬ 
parison  with  SHG  have  been  ma de, 21,22  but  only  when  absorption  and  dispersion 

are  negligible  (i.e.,  the  same  conditions  for  Klei nman's  symmetry17).  These 

c  cur 

comparisons  should  not  be  made  directly  between  r^  and  .  Instead,  they 
should  be  done  between  the  electro-optic  susceptibility,19 


2  2  S 
-n  n  r.. 
eo  o  e  41 
X41 


-1.5  x  10‘12  m/V 


(14) 


SHG 

and  the  second  harmonic  susceptibility,  d^  .  Since  a  non-zero  SHG  coef¬ 
ficient  for  Te02  violates  Klei nman's  symmetry,  the  necessary  assumptions  about 
absorption  and  dispersion  must  be  inappropri ate  and,  hence,  the  relationship 
between  the  electro-optic  and  SHG  coefficients  for  TeC^  is  not  straight¬ 
forward. 

Qualitatively,  the  SHG  coefficient,  while  being  non-zero,  is  still 
recognized  as  being  small  in  magnitude  when  compared  to  other  materials.  As  a 
consequence,  the  electronic  contribution  to  the  electro-optic  effect  is 
assumed  to  be  approximately  zero.  Then,  since  the  measured  electro-optic 
coefficient  is  also  small,  the  conclusion  is  that  the  expected  soft-mode 
enhancement  of  the  ionic  contributions  to  the  electro-optic  susceptibility 
does  not  occur. 
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Crystal  orientation  for  Te02  electro-optic  sample. 

Optical  configuration  of  electro-optic  measurement. 

Interference  fringe  patterns  in  Te02  for  (a)  -4.9,  (b)  0.0  and  (c) 
+4.9  kV. 
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